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The  Navy’s  mine  impact  burial  prediction  model  creates  a  time  history  of  a  cylindrical  or 
a  noncylindrical  mine  as  it  falls  through  air,  water,  and  sediment.  The  output  of  the  model 
is  the  predicted  mine  trajectory  in  air  and  water  columns,  burial  depth/orientation  in 
sediment,  as  well  as  height,  area,  and  volume  protruding.  Model  inputs  consist  of  param¬ 
eters  of  environment,  mine  characteristics,  and  initial  release.  This  paper  reviews  near 
three  decades  ’  effort  on  model  development  from  one  to  three  dimensions:  (1)  one¬ 
dimensional  models  predict  the  vertical  position  of  the  mine’s  center  of  mass  (COM)  with 
the  assumption  of  constant  falling  angle,  (2)  two-dimensional  models  predict  the  COM 
position  in  the  (x,z)  plane  and  the  rotation  around  the  y-axis,  and  (3)  three-dimensional 
models  predict  the  COM  position  in  the  (x,y,z)  space  and  the  rotation  around  the  x-,  y-, 
and  z-axes.  These  models  are  verified  using  the  data  collected  from  mine  impact  burial 
experiments.  The  one-dimensional  model  only  solves  one  momentum  equation  (in  the 
z-direction).  It  cannot  predict  the  mine  trajectory  and  burial  depth  well.  The  two- 
dimensional  model  restricts  the  mine  motion  in  the  (x,z)  plane  (which  requires  motion¬ 
less  for  the  environmental  fluids)  and  uses  incorrect  drag  coefficients  and  inaccurate 
sediment  dynamics.  The  prediction  errors  are  large  in  the  mine  trajectory  and  burial 
depth  prediction  (six  to  ten  times  larger  than  the  observed  depth  in  sand  bottom  of  the 
Monterey  Bay).  The  three-dimensional  model  predicts  the  trajectory  and  burial  depth 
relatively  well  for  cylindrical,  near-cylindrical  mines,  and  operational  mines  such  as 
Manta  and  Rockan  mines.  [DOI:  10.1115/1.3013823] 
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1  Introduction 

In  mine  hunting,  success  often  hinges  on  knowing  as  much  as 
possible  about  the  mines  that  have  been  placed  and  the  environ¬ 
ment  that  has  affected  the  placement.  Since  bottom  mines  cannot 
be  searched  visually,  and  are  often  difficult  to  locate  with  conven¬ 
tional  sonar,  an  estimate  of  area  and  height  of  the  mine  protruding 
from  the  sediment,  or  the  burial  depth,  is  crucial  for  the  planning 
and  execution  of  mine  clearance  operations.  Determining  the 
likely  mine  burial  depth  requires  numerical  models  of  the  burial 
process  and  knowledge  of  the  environment,  including  sediment 
properties,  waves,  tides,  and  water  depth. 

Sea  deployed  mines  currently  used  by  the  United  States  and 
other  nations  fall  into  three  general  categories:  bottom,  moored, 
and  drifting  mines.  Bottom  mines  rest  on  the  ocean  floor  and  are 
generally  deployed  in  littoral  regions.  Common  placements  for 
bottom  mines  include  shipping  channels,  harbors,  anchorages,  riv¬ 
ers,  and  estuaries.  Bottom  mines  are  deployed  in  one  of  the  three 
ways:  aircraft,  surface  ship,  and  submarine.  Mine  impact  burial 
models  have  been  developed  to  predict  mine’s  motion  in  air  and 
water  and  to  determine  the  burial  depth  in  sediment  when  the 
mine  comes  to  rest  in  the  sediment. 

One-dimensional  impact  burial  prediction  model  (IBPM)  was 
developed  by  Arnone  and  Bowen  [1]  to  predict  the  vertical  posi¬ 
tion  of  the  cylindrical  mine’s  center  of  mass  (COM)  as  it  falls 
through  air,  water,  and  sediment.  The  burial  depth  of  the  mine  in 
marine  sediment  is  then  calculated  from  the  mine’s  velocity  and 
the  sediment  characteristics.  IBPM  only  solves  the  vertical  mo¬ 
mentum  equation  with  the  assumption  of  an  unchanged  orienta¬ 
tion  in  the  fluid. 
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Satkowiak  [2,3]  advanced  Arnone  and  Bowen’s  [1]  pioneering 
work  through  correcting  reference  flow  for  drag  and  added-mass 
calculations;  reworking  equations  for  sediment-cavity  regime, 
drag  due  to  the  cylindrical  and  rounded  noses,  and  resistant  forces 
in  semisolid  sediment;  and  including  water  temperature  effect  on 
the  water  viscosity.  The  major  weakness  of  the  one-dimensional 
model  is  the  mine’s  orientation  (or  the  falling  angle)  assumed 
constant  as  it  falls  through  the  fluid. 

Two-dimensional  models  were  developed  first  by  Hurst  to  over¬ 
come  the  major  weakness  of  IBPM  (constant  falling  angles)  [4]. 
The  models,  written  in  basic  (IMPACT25)  and  in  matlab  (IM- 
PACT28),  contain  two  momentum  equations  (in  the  x-  and 
z-directions)  and  one  moment-of-momentum  equation  (in  the 
y-direction),  include  Mulhearn’s  [5]  formulation  for  sediment 
bearing  strength,  and  use  multilayered  sediments.  They  predict 
mine  COM  position  in  the  (x,z)  plane  and  the  rotation  (i.e.,  mine 
orientation)  around  the  y-axis. 

Although  the  two-dimensional  models  advance  our  knowledge 
on  the  mine  movement  by  including  its  rotation  around  the  y-axis, 
it  is  very  difficult  to  include  the  motion  of  the  fluid.  This  is  be¬ 
cause  it  is  hard  to  assume  the  fluid  (air,  water,  and  sediment) 
movement  strictly  in  the  (v,z)  plane.  Any  fluid  motion  in  the 
y-direction  induces  drag  force  and  in  turn  causes  the  mine  move¬ 
ment  (in  the  y-direction),  which  breaks  the  two-dimensional  sce¬ 
nario  [6,7].  In  fact,  it  is  impossible  to  lay  a  mine  in  the  same 
vertical  plane  of  the  fluid  velocity.  Sensitivity  studies  on 
IMPACT25/28  burial  depth  show  the  insensitivity  of  the  mine 
releasing  height  and  the  water  temperature  [8].  A  mine  drop  ex¬ 
periment  at  Monterey  Bay  [9-11]  shows  that  IMPACT25/28  over¬ 
predicts  the  burial  depth. 

Three-dimensional  model  (IMPACT35)  has  been  developed 
with  the  support  from  the  Office  of  Naval  Research  through  Mine 
Burial  Prediction  Program  [12-15].  The  model  contains  three  mo¬ 
mentum  equations  and  three  moment-of-momentum  equations.  It 
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Fig.  1  M-coordinate  with  the  COM  as  the  origin  of  X  and 
(imJm)  as  the  two  axes.  Here,  \  is  the  distance  between  the 
COV  (B)  and  Com  (X),  and  ( L,R )  are  the  cylinder’s  length  and 
radius  (after  Chu  et  al.  [13]). 


predicts  the  mine’s  COM  position  in  the  (x,y,z)  space  and  the 
rotation  (i.e.,  mine’s  orientation)  around  the  three  axes.  Several 
mine  drop  experiments  conducted  at  the  Naval  Postgraduate 
School  (NPS),  Naval  Undersea  Warfare  Center  (NUWC)- 
Carderock,  and  Baltic  Sea  (by  the  German  Navy)  were  used  to 
evaluate  the  two-  and  three-dimensional  models.  The  results  show 
great  improvement  of  the  three-dimensional  modeling. 

The  one-,  two-,  and  three-dimensional  models  are  reviewed  in 
this  paper.  Basic  physics,  formulation,  strength,  and  weakness  of 
each  model  are  presented.  The  purpose  is  to  provide  a  overview  of 
more  than  two  decades’  effort  on  predicting  the  mine  (cylindrical 
and  operational  mine  shape)  movement  in  the  air,  water,  and  sedi¬ 
ment  column. 

2  Mine  Location  and  Orientation 

Consider  an  axially  symmetric  cylinder  with  the  COM  X  (or 
called  gravity  center  (GC)  in  literatures)  and  the  center  of  volume 
(COV)  B  on  the  main  axis  (Fig.  1).  Let  (L,R,x)  represent  the 
cylinder  length,  the  radius,  and  the  distance  between  the  two 
points  (X,B).  The  positive  ^-values  refer  to  the  nose-down  case, 
i.e.,  point  X  is  lower  than  point  B.  Let  FE(0, i , j , k)  be  the  earth- 
fixed  coordinate  (F-coordinate)  with  the  origin  O,  and  three  axes: 
x-,  y-axes  (horizontal)  with  the  unit  vectors  (i,j)  and  z-axis  (ver¬ 
tical)  with  the  unit  vector  k  (upward  positive).  The  position  of  the 
cylinder  is  represented  by  the  position  of  the  COM, 


X  =  xi  +  yj  +  zk  (1) 

which  is  a  translation  of  the  cylinder.  The  translation  velocity  is 
given  by 

dX 

—  =  V,  \  =  (u9v,w)  (2) 

dt 

The  orientation  of  the  cylinder’s  main  axis  (pointing  down¬ 
ward)  is  given  by  iM.  The  angle  between  \M  and  k  is  denoted  by 
02  +  7i/  2  (Fig.  2).  The  angle  02  is  the  mine  falling  angle.  In  the 
one-  and  two-dimensional  modelings,  only  the  F-coordinate  sys¬ 
tem  is  used.  In  three-dimensional  modeling,  two  extra-coordinate 
systems  (main-axis  following  and  force  following  coordinates)  are 
also  used. 

3  Triple  Coordinate  Systems 

Three  coordinate  systems  are  used  in  mine  impact  burial  pre¬ 
diction  modeling:  earth-fixed  coordinate  (F-coordinate),  main- 
axis  following  coordinate  (M-coordinate),  and  force  following  co¬ 
ordinate  (F-coordinate)  systems.  All  the  three  coordinate  systems 
are  three  dimensional,  orthogonal,  and  right  handed  [16].  Projec¬ 
tion  of  the  mine’s  main-axis  vector  \M  onto  the  (v,y)  plane  creates 
angle  (03)  between  the  projection  and  the  v-axis  (Fig.  2).  The 
M-coordinate  is  represented  by  FM(X,iM,jM,kM)  with  the  origin 
X  (i.e.,  the  COM  location),  unit  vectors  (iM,jM,kM),  and  coordi¬ 
nates  ( xM,yM,ZM )•  The  unit  vectors  of  the  M-coordinate  system 
are  given  by 

Sm  -  k  X  iM,  kM  =  iM  X  jM  (3) 

The  M-coordinate  system  is  solely  determined  by  orientation  of 
the  cylinder’s  main  axis  \M. 

The  F-coordinate  is  represented  by  FF(X,iFJF,kF)  with  the 
origin  X,  unit  vectors  (iF,jF,kF),  and  coordinates  (xF,yF,zE)-  Let 
Vw  be  the  fluid  velocity.  The  water-to-mine  velocity  is  represented 
by 

Vr  =  vw-V  (4) 

which  can  be  decomposed  into  two  parts 

Vr  =  V1  +  V2,  V,  =  (Vr-i,.)i,,  \2  =  \r-(\r-iF)iF  (5) 

where  \i  is  the  component  paralleling  to  the  cylinder’s  main  axis 
(i.e.,  along  iM),  and  V2  is  the  component  perpendicular  to  the 


Fig.  2  Mine’s  COM  position  ( x,z )  and  orientation  02  (after  Chu  et  al.  [16]) 
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cylinder’s  main-axial  direction.  The  unit  vectors  for  the 
F-coordinate  are  defined  by  (column  vectors) 

=  jF-V2/|V2|,  —  i/r  X  (6) 

Transforms  among  the  three  coordinate  systems  can  be  found  in 
Ref.  [16]. 

4  One-Dimensional  Modeling 

One-dimensional  models  assume  that  the  cylinder  is  not  rotat¬ 
ing  (i.e.,  constant  orientation)  about  any  axis,  nor  a  net  fluid  dy¬ 
namic  lift.  Consequently,  it  can  only  be  applied  to  stable  motion 
of  the  body  along  one  of  its  major  axes  (horizontal  or  vertical). 
The  models  predict  the  COM  location  (i.e.,  z)  using  the 
Kirchhoff-Kelvin  theory.  Only  the  F-coordinate  system  is  used 
for  the  one-dimensional  modeling. 

4.1  Kirchhoff-Kelvin  Theory.  In  the  later  half  of  the  19th 
century,  Kirchhoff  and  Kelvin  showed  that  the  motion  of  a  solid 
moving  through  an  ideal  fluid  could  be  represented  by  a  compact 
system  of  equations  describing  the  coupled  fluid-body  dynamics 
[17,18].  Publications  about  the  implications  and  extensions  of  this 
theory  (cited  as  the  Kirchhoff-Kelvin  equations  hereafter)  were 
frequent  from  the  1870s  until  the  1900s,  when  experimental  aero¬ 
dynamics  showed  the  limitations  of  the  idealized  models  in  rep¬ 
resenting  the  coupled  dynamics  of  bodies  and  complex  turbulent 
boundary  layer  regimes.  Interest  in  the  subject  also  waned  due  to 
the  daunting  nonlinearity  of  the  governing  equations.  Recently 
there  has  been  a  resurgence  of  interest  in  the  Kirchhoff-Kelvin 
theory,  since  the  widespread  availability  of  highly  capable  com¬ 
puters  has  allowed  numerical  analysis  of  problems  involving  non¬ 
linear  dynamics  [19-21].  Even  when  analytical  closure  is  not 
practical,  the  general  results  of  the  Kirchhoff-Kelvin  equations 
aid  heuristic  explanations  of  the  observed  dynamics  of  more  com¬ 
plex  shapes,  such  as  coins  (very  short  cylinders)  sinking  in  water 
[22]  and  tumbling  cards  [23].  In  Ref.  [24],  excellent  introductory 
can  be  found  to  the  Kirchhoff-Kelvin  equations  and  implications. 

In  general,  the  Kirchhoff-Kelvin  equations  predict  the  coupled 
dynamic  response  of  arbitrary  solid  bodies  to  various  forces  and 
torques  within  inviscid  incompressible  flow.  For  high  Reynolds 
number  regimes  (R q=UD/  v>  104,  where  U  is  the  velocity,  D  is 
either  the  length  or  diameter  of  the  cylinder  depending  on  release 
orientation,  and  v  is  the  kinematic  viscosity),  the  direct  effect  of 
viscosity  is  small  such  that  coherent  structures  of  the  turbulent 
boundary  layer  diminish.  The  typical  value  of  Re  for  mine  falling 
through  the  water  column  is  around  105  [15].  This  condition  al¬ 
lows  the  generalized  dynamics  of  falling  bodies  to  be  character¬ 
ized  in  terms  of  a  simplified  form  of  the  Kirchhoff-Kelvin  theory 
where  buoyancy  force  and  turbulent  drag  are  balanced  by  the 
inertia  of  the  cylinder  and  displaced  water.  For  a  freely  sinking 
cylinder  of  diameter  d  and  length  /,  the  combined  effect  of  buoy¬ 
ancy  and  gravity  per  unit  volume  is 

B  =  —  (p  —  pw)g  (7) 

where  g  is  the  gravitational  acceleration,  and  ( p-pw )  is  the  den¬ 
sity  difference  of  the  cylinder  and  water.  The  flow  over  the  cylin¬ 
der  surface  yields  pressure  distributions  that  result  in  a  net  drag 
force  (Fd)  on  the  body,  empirically  represented  as  proportional  to 
the  square  of  the  body’s  speed  through  water, 

Fd~\PwCdU2  A  (8) 

where  Cd  is  the  drag  coefficient,  U  is  the  mine’s  falling  speed,  and 
A  is  an  effective  cross-sectional  area  normal  to  the  flow.  Under 
the  assumption  that  the  torques  exerted  by  wake  pressure  fluctua¬ 
tions  are  insignificant,  the  phenomenological  representation  of  the 
total  drag  force  on  the  cylinder  can  be  used  in  an  approximation 
for  the  falling  motion  where  the  rate  of  change  in  kinetic  energy  is 
parametrized  by  the  rate  of  change  of  the  cylinder  momentum  and 
the  acceleration  of  the  “virtual  mass”  of  the  displaced  water 


(which  depends  upon  the  geometric  characteristics  of  the  body). 
Thus,  the  force  balance  on  the  cylinder  is  given  by 

_  PWQA  /  dz  \  2  _  ip  -  pj  +  F_ 
dt2  2pXl{\  +  f)\ dt }  p  g  pH 

where  IT  is  the  volume  of  the  mine,  /  is  the  effective  added-mass 
factor  or  the  virtual  mass  coefficient  of  the  system  due  to  the 
acceleration  of  water  around  the  moving  body,  and  F  is  the  addi¬ 
tional  sediment  force,  which  is  zero  in  air  and  water.  The  left-hand 
side  of  Eq.  (9)  shows  the  vertical  acceleration  of  the  mine  per  unit 
mass.  The  right-hand  side  of  Eq.  (9)  shows  the  external  forces  per 
unit  mass  exerted  on  the  mine,  with  the  first  term  from  the  drag 
force  (Fd),  the  second  term  from  the  buoyancy  force  (F),  and  the 
third  term  from  the  sediment  resistance. 

Full  implementation  of  the  theory  is  not  the  focus  of  this  work; 
however,  stability  analysis  of  particular  variants  of  the  Kirchhoff- 
Kelvin  equation  provides  useful  insights.  We  can  characterize  the 
generalized  dynamics  governing  the  motions  of  mines  falling 
through  water  under  the  following  hypothetical  situations:  vertical 
descent  with  the  axis  of  the  cylinder  aligned  parallel  to  the  flow  or 
vertical  descent  with  the  axis  of  the  cylinder  aligned  perpendicular 
to  the  flow.  For  the  simplest  case  where  the  axis  of  the  cylinder  is 
aligned  with  the  fall  direction,  body  motions  will  be  stable  unless 
they  are  significantly  disturbed  in  a  direction  normal  to  the  mo¬ 
tion.  This  means  that  a  vertically  oriented  cylinder  dropped  freely 
will  maintain  this  orientation  until  external  forces,  such  as  vari¬ 
able  turbulent  drag,  cause  the  motion  to  become  unstable  resulting 
in  changes  in  orientation.  With  the  axis  of  the  cylinder  oriented 
horizontally,  or  normal  to  the  trajectory,  the  motion  is  also  stable. 
An  alternative  scenario  where  the  axis  of  the  cylinder  is  not  ini¬ 
tially  aligned  either  parallel  or  perpendicular  to  the  flow  makes  the 
application  of  the  above  equations  difficult  but  has  been  addressed 
in  Ref.  [25],  who  established  that  with  time,  the  falling  body  will 
assume  a  horizontal  orientation  where  the  torque  exerted  on  the 
body  by  the  water  turns  the  main  axis  of  the  cylinder,  so  that  it  is 
normal  to  the  relative  flow  (“broadside”).  This  configuration  is 
stable  when  the  exposed  area  is  normal  to  the  flow  and  is  signifi¬ 
cantly  greater  than  the  cross-sectional  area  of  the  “nose-on”  atti¬ 
tude.  Slight  variations  of  this  stable  mode  also  exist  including 
regular  oscillations,  glide  tumble  motions,  and  helical  motions. 

4.2  IBPM.  The  one-dimensional  model  is  the  first  generation 
of  the  Navy’s  mine  impact  burial  model  [1-3],  called  the  IBPM. 
This  model  is  used  to  predict  the  vertical  location  and  orientation 
of  a  cylindrical  mine  falling  through  air,  water,  and  sediment 
through  solving  the  Kirchhoff-Kelvin  equation  (9).  The  model 
consists  of  four  major  components:  (a)  steady  falling  attitude  (02), 
(b)  drag  computation,  (c)  cavity  regimes,  and  (d)  sediment  forces. 


4.2.1  Steady  Falling  Angle  in  Single  Media.  IBPM  does  not 
use  the  moment-of-momentum  equation  for  predicting  mine  ori¬ 
entation.  Instead,  the  tangent  of  the  falling  angle  (or  called  atti¬ 
tude)  (02)  is  assumed  as  the  ratio  of  moments  due  to  the  axial  (Fa) 
and  cross  (Fc)  forces  [1]: 


tan  02  - 


moment(Fj 

moment(Fc) 


(10a) 


where  the  axial  and  cross  forces  are  calculated  by 

Fa  =  g(p-fiPw  sin2  02).  Fc  =  g(p-f2pw  COS2  ip2)  (10/?) 

Here,  (/j  ,/2)  are  the  effective  added-mass  factors  in  the  axial  and 
cross  mine  directions, 


fi  = 


-  tfo 


fl  - 


A) 
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The  two  parameters  (cr0,/30)  depend  only  on  the  aspect  ratio  (£), 
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1  In 


2+8 

2^8, 


(12) 


The  effective  added-mass  factor  is  computed  by 

f=h  sin  <f>2  +/i  cos  i/j2  (13) 

For  an  attitude  02  of  0  deg  (vertical),  /  equals  /j.  For  an  attitude 
02  of  90  deg  (horizontal),  /  equals  f2.  In  IB  PM,  the  mine’s  atti¬ 
tude  (02)  is  obtained  by  solving  Eqs.  (10a)  and  (10&)  using  the 
iteration  method.  Since  the  expression  (10&)  uses  the  fluid  density 
(pw),  the  mine’s  attitude  (02)  does  not  change  in  a  single  fluid  but 
does  change  when  it  passes  through  the  interface  of  two  fluids. 
The  effective  cross-sectional  area  (or  sometimes  called  projected 
area)  A  depends  on  the  attitude  02, 

A  =  2LR  sin  ip2  +  ttR 2  cos  ip2  (14) 

which  is  used  for  the  drag  calculation  (see  the  next  subsection). 
For  an  attitude  of  0  deg  (vertical),  A  equals  ttR2.  For  an  attitude 
of  90  deg  (horizontal),  A  is  the  product  of  the  length  (L)  and 
diameter  (2 R).  It  is  noted  that  the  physical  base  for  “constant  02  in 
a  single  fluid”  is  weak.  This  is  because  the  mine’s  attitude  (02)  is 
determined  by  solving  the  moment-of-momentum  equation  (see 
Secs.  5.1  and  6.1)  rather  than  solving  Eqs.  (10a)  and  (10&). 


4.2.2  Drag  Coefficient.  Drag  is  caused  by  hydrodynamic 
forces  acting  on  the  falling  cylinder  in  the  axial  and  cross  direc¬ 
tions,  and  therefore  it  depends  on  the  attitude  (02)  of  the  cylinder. 
The  drag  becomes  somewhat  more  complicated  because  the  fea¬ 
ture  of  the  flow  across  the  cylinder  changes  based  on  the  magni¬ 
tude  of  the  Reynolds  number, 


Re  = 


2  UR 


(15) 


where  v  is  the  kinematic  viscosity  of  the  fluid.  When  the  Reynolds 
number  is  small  (<1),  skin  friction  dominates  but  as  the  Reynolds 
number  increases,  a  more  laminar  flow  occurs  and  the  pressure 
drag  dominates.  At  a  certain  and  critical  Reynolds  number  (104), 
flow  becomes  turbulent  and  drag  abruptly  decreases.  Reynolds 
number  is  thus  dependent  on  velocity  of  the  cylinder  and  must  be 
constantly  recalculated  as  the  cylinder  falls  through  each  medium. 

Dependence  of  drag  on  (Re,  if/2)  is  shown  in  the  axial  and  cross 
drag  coefficients  ( Cda,Cdc )  [1] 

Cdc  =i.lrjCda3  sin2  t//2,  Cda=\.\CdifA-  (16) 


where 


J 


cd°°~* 


0.84864  +  5.81939/Re 
(0.5833  X  10"4)Re  + 0.61677 
1.2 

1.19381  +  (0.65828  X  10"6)Re-  10“u 
8  X  10-7  Re -0.15 


Re2 


if  Re  <  4  X  103 
if  4  X  103  <  Re  <  104 
if  104  <  Re  <  8  X  104 
if  8  X  104  <  Re  <  5  X  105 
if  Re  &  5  X  105 


(17) 


r 


is  the  cross  flow  drag  coefficient  for  the  cylinder  with  infinite 
length  (Cd oo),  and 


Q*  -  0.33£+ 


3.984 


VRe  _ 


-  +  ,51/2 
S 


(18) 


is  the  axial  drag  coefficient  with  no  surface  imperfections.  The 
coefficient  of  1.1  in  Eq.  (16)  is  used  to  account  for  imperfections 
of  the  cylinder  with  the  correction  factor  77  given  by 

1  (A2  ,  ^ 

7]  =  0.52238  +  0.021 19-  -  0.000481  -  I  (19) 


The  total  drag  coefficient  (Cd)  for  computing  the  drag  force  on  a 
mine  moving  through  the  fluid  is  calculated  by 

Q  -  Qa  +  Cdc  Sin2  02  (20) 


trend  of  this  ratio  is  known.  For  example,  the  void  ratio  equals  “0” 
for  the  cylinder  in  the  air  and  equals  “1”  for  the  cylinder  totally 
wetted  in  the  water. 

When  the  pressure  generated  by  the  dynamic  loading  of  the 
cylinder  impacting  on  the  water  is  balanced  by  the  hydrostatic 
pressure  of  the  water.  The  air-water  cavity  is  assumed  to  collapse 
and  the  cylinder  is  fully  wetted.  The  dynamic  pressure  generated 
by  the  cylinder  ( pD )  is  the  sum  of  the  atmospheric  pressure  pa 
(=100  kPa)  at  the  water  surface  and  dynamic  pressure  due  to  the 
cylinder’s  falling  velocity  U, 

PD(t)=Pa+\Pw\U{t)Y  (21  a) 


The  hydrostatic  pressure  ( ps )  of  water  is  the  sum  of  the  atmo¬ 
spheric  pressure  ( pa )  and  weight  of  the  water  column, 


4.2.3  Air-Water  Cavity.  Upon  impacting  the  water,  the  cylin¬ 
drical  mine  enters  the  air-water  cavity.  The  properties  of  the  fluid 
in  the  water  cavity  are  combination  of  air  and  water  fluid  proper¬ 
ties  and  are  continuously  changing  with  time.  It  is  extremely  dif¬ 
ficult  to  accurately  predict  the  forces  acting  on  a  body  in  a  fluid  of 
changing  properties  [1,26].  The  changing  properties  around  the 
cylinder  through  the  cavity  are  the  fluid  density  and  kinematic 
viscosity.  For  example,  the  fluid  density  changes  from  air  ( pa 
=  1.29  kgm-3)  to  sea  water  (pw=1025  kgm-3). 

Within  the  air-water  cavity  regime,  a  percentage  of  each  of 
these  densities  is  used  in  determination  of  the  resulting  average 
density.  This  is  represented  by  the  ratio  (called  the  void  ratio)  of 
the  volume  of  water  in  the  cavity  to  the  total  cavity  volume. 
Although  it  is  difficult  to  actually  determine  these  volumes,  the 


Ps(t)=Pa~  Pwgzit)  (21  b) 


In  the  IB  PM,  the  square  root  of  the  ratio  between  pD  and  ps  is 
defined  as  the  void  ratio  n0 


n0(t )  = 


(22) 


which  is  taken  as  a  weight  of  water  in  the  air-water  cavity.  The 
cavity  density  (p^av)  and  viscosity  (z/^av)  are  calculated  by 


PT=PaV  -«o)  +  Pw«0> 


ncav  cav 

r'w  w 


(1  -«o)  »0 

Pava  Pwvw 


(23) 
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The  dynamic  effect  of  cavitation  on  the  body  is  through  the 
change  of  the  Reynolds  number  (in  turn  the  change  of  the  drag 
coefficient).  The  Reynolds  number  for  the  air- water  cavity  is  cal¬ 
culated  by 


Re 


cav 

w 


2  UR 


Substitution  of  Re^av  into  Eqs.  (16)-(20)  leads  to  the  drag  coeffi¬ 
cient  for  the  air-water  cavity. 


4.2.4  Water-Sediment  Cavity.  Upon  impacting  the  sediment, 
the  cylindrical  mine  creates  a  cavity  in  which  the  fluid  properties 
of  water  and  sediment  are  interacting.  Similar  to  the  air-water 
cavity  regime,  the  properties  of  the  water- sediment  cavity  are 
combination  of  water  and  sediment  properties  and  are  continu¬ 
ously  changing  with  time.  In  the  water- sediment  cavity  regime, 
the  fluid  density  changes  from  water  (pw=  1025  kgm-3)  to  sedi¬ 
ment  ps.  Within  the  water- sediment  cavity,  the  void  ratio  nx 


is  taken  as  the  weight  of  the  sediment  in  water- sediment  cavity. 
Here, 


Pd(0  =  Pa  +  Pwgh  +  \ps\.U{t)f,  psit )  =  Pa  +  pwgh  -  psg[z(t )  +  h] 

(25) 


are  the  dynamic  and  static  pressures,  and  h  is  the  water  depth.  The 
water- sediment  cavity  density  (pjav)  and  viscosity  (i^av)  are  cal¬ 
culated  by 


pT'  =  pJ1  -«i)  +  p.v«b 


i  (i-«i)  »i 

pT^T~  w*  +  Ps** 


(26) 


The  kinematic  viscosity  of  the  sediment  is  determined  by 


v*  =  vw  +  ■ 


psdU/dz 


(27) 


where  S  is  the  shear  strength  of  the  sediment.  Substitution  of  vs 
into  Eq.  (15)  leads  to  a  new  Reynolds  number;  and  then  use  of 
Eqs.  (16)-(20)  leads  to  the  drag  coefficient  for  the  water-sediment 
cavity. 


4.2.5  Sediment  Forces.  Forces  exerting  on  mine  in  sediment 
are  different  from  that  in  air  and  water,  especially  the  sediment 
resistant  force  (F),  which  is  calculated  in  IB  PM  using  two  differ¬ 
ent  methods.  The  first  method  [1]  is  velocity  independent.  It  de¬ 
composes  the  sediment  resistant  force  into  compressive  force  (Fc) 
and  shearing  force  (F5), 

F  =  FC  +  FS  (28  a) 

where  (FC,FS)  are  proportional  to  the  shear  strength  (5)  with  the 
proportionalities  of  the  cylindrical  reference  flow  areas  (Af,As)  in 
the  front  (twice)  and  in  the  side  (once), 

Fc  =  2SAf,  Fs  =  SAs  (28  b) 

The  second  method  is  velocity  dependent  [3].  The  sediment  resis¬ 
tant  force  is  calculated  by 

F  =  EsS(NbAf+  SafAs/St)  (29) 

where  Es  is  the  sediment  strain  rate  (i.e.,  ratio  of  the  sediment 
strength  for  a  given  velocity  of  the  cylinder),  Nb  (=10)  is  the 
bearing  capacity  factor,  5af  is  the  side  adhesive  factor  (0.3  for  the 
cavity  and  1.0  for  sediment),  and  St  (=3)  is  the  sediment  sensitiv¬ 
ity  coefficient. 

Thus,  the  sediment  shear  strength  and  density  data  are  the  in¬ 
puts  to  the  Kirchhoff-Kelvin  equation  (9)  for  mine  movement  in 
the  sediment  since  they  are  needed  for  determination  of  the  sedi¬ 


ment’s  compressive  and  shearing  forces,  buoyancy  force,  added 
mass,  drag  force,  and  kinematic  viscosity.  Since  the  properties  of 
sediment  change  with  depth,  the  sediment  density  and  shear 
strength  profiles  with  5  cm  vertical  resolution  are  used  in  the 
IBPM. 


4.3  Sensitivity  Studies.  For  a  single  medium,  after  the  atti¬ 
tude  (02)  is  determined  from  Eq.  (10a)  and  (10/?),  the  effective 
cross-sectional  area  A  and  the  effective  added-mass  factor  /  can 
be  calculated  using  Eqs.  (11)  and  (12),  and  then  the  total  drag 
coefficient  Cd  can  be  computed  using  Eq.  (20).  For  sediment,  the 
compressive  and  shear  stress  forces  (FC,FS)  are  computed  from 
the  density  and  shear  strength.  With  the  known  parameters  (A,  /, 
Cd,  Fc,  and  Fs),  the  Kirchhoff-Kelvin  equation  (9)  can  be  solved 
for  mine  movement.  For  mine  penetration  into  the  air-water  and 
water- sediment  interfaces,  the  total  drag  coefficient  Ccdy  is  calcu¬ 
lated  from  the  cavity  density  (p^av  or  pjav)  and  viscosity  (jF*v  or 

*n- 

Arnone  and  Bowen  [1]  conducted  a  model  sensitivity  study  on 
the  mine  air  weight  (or  wet  weight),  length,  and  radius.  The  model 
(9)  was  integrated  with  various  combinations  of  these  parameters 
and  three  different  attitudes:  02-0  deg,  (horizontal),  45  deg,  and 
90  deg  (vertical).  The  water  impact  velocity  has  the  following 
features,  (a)  For  mine  movement  in  the  air,  as  one  would  expect, 
the  more  streamlined  the  falling  attitude  (i.e.,  0 2-90  deg  verti¬ 
cal),  the  higher  velocity  the  cylindrical  mine  falls,  (b)  The  water 
impact  velocity  varies  drastically  with  the  attitude  for  light  mines. 
Such  an  effect  reduces  as  the  mass  increases,  (c)  The  water  impact 
velocity  is  not  sensitive  to  mine  length  (L)  and  radius  ( R ). 

The  mine  falling  velocity  in  the  water  has  the  following  char¬ 
acteristics.  (a)  It  varies  drastically  with  attitude  and  wet  weight. 
For  the  same  wet  weight,  the  mine  falls  faster  for  02-90  deg 
(vertical)  than  for  02  =  0  deg  (horizontal).  For  the  same  02,  the 
mine  falls  faster  for  heavier  wet  weight,  (b)  It  is  not  sensitive  to 
attitude  for  short  mines  and  is  very  sensitive  to  attitude  for  long 
mines,  (c)  It  is  not  sensitive  to  mine  length  for  horizontal  release 
(02  =  0  deg)  and  is  very  sensitive  to  mine  length  for  vertical  re¬ 
lease  (02  =  90  deg),  (d)  It  is  not  sensitive  to  attitude  for  small  R 
and  is  very  sensitive  to  attitude  for  large  R.  (e)  It  is  not  sensitive 
to  attitude  for  small  R  and  is  very  sensitive  to  attitude  for  large  R. 

4.4  Strength  and  Weakness.  IBPM  has  capability  to  simulate 
mine  falling  velocity  in  the  air  and  water  columns.  As  mentioned 
in  the  previous  subsection,  the  model  provides  the  following  use¬ 
ful  results:  (a)  The  mine  has  higher  falling  velocity  for  vertical 
release  (02= 90  deg)  than  for  the  horizontal  release  (02  =  0  deg), 
(b)  the  water  impact  velocity  varies  drastically  with  the  attitude 
for  light  mines  but  does  not  for  heavy  mines,  and  (c)  the  mine 
falling  velocity  in  water  column  is  sensitive  to  attitude,  wet 
weight,  mine  length,  and  mine  radius.  For  the  same  mine,  the 
falling  velocity  in  water  is  minimum  for  horizontal  orientation  and 
maximum  for  vertical  orientation. 

Weakness  of  the  IBPM  is  to  assume  the  constant  attitude  (02) 
during  the  mine  falling  through  a  single  medium  (air,  water,  or 
sediment)  (Fig.  3),  which  is  physically  unrealistic  [12].  In  fact,  the 
orientation  changes  for  any  solid  object  falling  through  fluid  (air 
and  water).  When  COV  is  not  colocated  with  COM,  the  buoyancy 
force  (B)  has  a  moment  of  momentum  exerted  on  mine  (Fig.  4) 
[27].  If  the  mine  slants  toward  left  (right),  the  moment  of  momen¬ 
tum  due  to  the  buoyancy  force  will  rotate  the  mine  clockwise 
(anticlockwise).  This  causes  the  spiral  motion  for  mine  falling 
through  fluid. 


5  Two-Dimensional  Modeling 

To  overcome  the  major  weakness  of  the  one-dimensional  model 
(i.e.,  constant  attitude  02),  Hurst  [4]  modified  the  IBPM  model 
through  allowing  the  cylinder  to  move  both  vertically  and  hori¬ 
zontally  (in  the  (x,z)  plane)  as  well  as  rotation  about  the  y-axis. 
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Without  Moment 


(32) 


Fig.  3  Mine’s  orientation  is  assumed  constant  by  the  one¬ 
dimensional  model  when  it  falls  through  a  single  fluid  (after 
Chu  [12]) 


d2ijj2  v 

J2-^  =  BX  00^2  + Ml 
at 

where  J2  is  the  moment  of  inertia  in  the  y-axis,  and  M\  is  called 
the  braking  torque  in  Ref.  [4]. 

5.2  IMPACT25/28.  The  two-dimensional  models,  usually 
called  IMPACT25  (written  in  basic)  and  IMPACT28  (written  in 
matlab),  are  the  second  generation  of  the  Navy’s  mine  impact 
burial  prediction  models  for  cylindrical  mines  and  were  based  on 
Eqs.  (30)-(32)  for  obtaining  (jc,z,^2)*  The  external  forcing  for 
IMPACT25/28  consists  of  drag  force  ( Fxh,Fzh )  and  braking  torque 
(Myh).  Since  the  mine  movement  is  restricted  in  the  (x,z)  plane,  it 
is  very  hard  to  include  the  motion  of  the  fluids  (air  or  water).  If 
the  fluid  has  velocity  in  the  y-direction,  the  mine’s  motion  cannot 
be  two  dimensional.  Thus,  in  IMPACT25  and  IMPACT28,  the 
fluid  (air  or  water)  is  assumed  motionless. 

5.2.1  Drag  Force.  The  drag  force  for  the  whole  cylinder  is 
calculated  in  two  directions:  along  the  cylinder’s  main  axis  (iM) 
and  across  the  cylinder  (kM), 

Fd  =  ~  Fdi^M  sin  >l>2  ~  cos  i fe)  (33) 

where  Fd  is  computed  using  the  drag  law 

Fd=\pwCd{u2  +  w2)h  (34) 


These  changes  mandated  the  use  of  more  complicated  dynamical 
system  than  the  IBPM.  Two  coordinate  systems  ( E -  and 
M-coordinates)  are  used  in  the  two-dimensional  models. 


5.1  Dynamical  System.  Let  the  cylinder  be  moving  in  the 
(x,z)  plane.  The  momentum  equations  in  the  (x,z)  directions  are 
given  by 


d2x  _  du  Fxh 
dt 2  dt  pll 


(30) 


d2z  _  dw  B  +  F\ 

dt 2  dt  8  p\\ 


(31) 


where  B  is  the  buoyancy  force,  and  (Fxh ,  Fzh)  are  the  components  of 
the  hydrodynamic  force  F^.  Since  the  cylinder  is  restricted  in  the 
(x,z)  plane,  the  only  possible  rotation  is  around  the  y-axis,  which 
is  described  by  the  attitude  ijj2.  The  moment-of-momentum  equa¬ 
tion  is  written  in  the  M-coordinate  system, 


Fig.  4  Change  of  mine  orientation  caused  by  momentum  due 
to  the  buoyancy  force  (after  Chu  [27]) 


The  drag  coefficient  Cd  is  computed  by 

Q-Q+Qn  (35) 


where  Cd  is  computed  in  the  same  way  as  in  the  one-dimensional 
model  (i.e.,  Eq.  (20)),  and  Cdn  is  the  drag  coefficient  for  the  nose. 


5.2.2  Braking  Torque.  The  braking  torque  Myh  in  Eq.  (32)  is 
calculated  by 

Myh  =  -lcdcRPwa>2L3Vc  (36) 


where 


Vc-u  cos  if/2  +  w  sin  if/2 


(37) 


is  the  cross-cylinder  velocity.  Equation  (36)  is  only  used  when  the 
mine  is  fully  immersed  in  a  single  fluid.  During  the  cavity  re¬ 
gimes,  a  different  calculation  is  applied.  Furthermore,  the  torque 
is  in  opposite  sign  to  rotation  and  thus  acts  as  the  brake  to  the 
rotation  of  mine.  For  a  single  fluid,  the  basic  equations  (30)-(32) 
are  integrated  in  association  with  Eqs.  (33)  and  (36). 


5.2.3  Sediment.  The  one-dimensional  models  (IPBM  and  its 
modifications)  treat  sediment  as  fluid  with  two  characterized  pa¬ 
rameters:  sediment  density  and  shear  strength.  Different  from  the 
one-dimensional  models,  the  two-dimensional  models  treat  the 
sediment  as  solid  that  undergoes  plastic  deformation  [4].  The  per¬ 
tinent  parameters  are  sediment  density  and  bearing  strength.  The 
bearing  strength  is  the  load  bearing  capacity  of  the  sediment  and 
defined  as  the  pressure  in  front  of  the  object  penetrating  the  sedi¬ 
ment.  It  is  related  to  the  shear  strength,  typically  larger  by  a  factor 
of  about  10.  Three  elements  contributing  to  the  sediment  resis¬ 
tance  on  the  penetration  of  a  falling  mine  are  included  in 
IMPACT25/28:  bearing  strength  of  the  sediment  (70%),  hydrody¬ 
namic  drag  (25%),  and  buoyancy  (5%)  [4]. 

As  the  mine  penetrates  the  sediment,  the  hydrodynamic  drag 
that  retards  the  penetration  is  calculated  by 


psA(u2  +  w2)  (  h 


ClL  +  Cl 


(38) 


where  C\  and  C2  are  the  drag  coefficients  for  low  Reynolds  num¬ 
ber  and  h  is  the  depth  in  sediment.  The  buoyancy  force  comes  into 
play  when  the  mine  impacts  the  sediment.  A  crater  is  formed  and 
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the  force  that  is  required  to  create  and  enlarge  the  crater  is  given 
by  Eq.  (4) 

B  =  TTgh(pw  -  p) A  (39) 

Two  assumptions  are  made  regarding  the  mine.  First,  the  buoy¬ 
ancy  force  acts  uniformly  on  the  portion  of  the  mine  in  contact 
with  the  sediment.  Second,  a  cavity  is  formed,  remains  at  the  aft 
end  of  the  mine,  and  leads  the  mine  surface  to  the  buoyancy  force. 
The  second  assumption  is  not  always  true.  For  some  cases,  the 
cavity  collapses  under  the  weight  of  the  sediment.  The  collapse 
limit  is  defined  and  approximated  as  the  depth  where  the  buoy¬ 
ancy  force  is  20%  of  the  bearing  force.  Besides,  the  drag,  buoy¬ 
ancy,  and  bearing  strength  of  sediment  all  depend  on  the  mine 
contact  area  A. 


5.2.4  Air-Water  Cavity.  As  mentioned  before,  when  the  mine 
penetrates  the  air-water  interface,  it  forms  a  cavity  behind  the 
mine  that  changes  the  drag  forces.  Throughout  the  cavity  regimes, 
nonsymmetric  forces  are  acting  on  the  mine  that  generate  torque 
and  affect  the  rate  of  rotation.  The  cavity  parametrization  is  dif¬ 
ferent  between  the  two-dimensional  models  (IMPACT25/28)  and 
one-dimensional  model  (IPBM).  The  cavity  formation  in  the  two- 
dimensional  models  is  controlled  by  the  cavitation  number  (/Vcav) 
defined  by 


N„ 


2(Pout-Pin) 

pw(u2  +  W 2) 


(40) 


where  (poxlt,pm)  are  the  hydrodynamic  pressure  outside  (i.e.,  wa¬ 
ter)  and  inside  the  cavity.  Pressure  difference  between  the  outside 
and  inside  of  the  cavity  increases  or  the  mine’s  velocity  decreases; 
the  cavitation  number  increases  until  the  eventual  collapse  of  the 
cavity.  The  drag  coefficient  of  the  cavity  is  the  function  of  the 
cavitation  number, 


Cd(Ncay)  =  Cd(l  +  bNcJ  (41) 

When  the  air- water  cavity  collapses,  7Vcav=0,  Cd( 0)  is  the  water 
drag  coefficient.  The  torque  in  the  air-water  cavity  can  be  com¬ 
puted  after  the  drag  coefficient  is  determined. 

5.3  Sensitivity  Study.  Sensitivity  studies  [8,28]  were  con¬ 
ducted  on  the  two-dimensional  model  (IMPACT25/28)  to  ascer¬ 
tain  which  parameters  the  model  is  most  sensitive  to  and  which 
can  be  eliminated  in  order  to  simplify  its  use.  The  model  was 
altered  to  allow  most  parameters  to  be  set  and  a  loop  run  for  one 
variable  at  a  time.  All  model  runs  were  made  with  preset  mine 
profile,  which  has  a  dry  weight  of  538  kg,  a  wet  weight  of  251  kg, 
and  a  uniform  diameter  of  0.475  m. 

5.3.1  Sensitivity  to  Release  Parameters.  Figure  5  demon¬ 
strates  the  sensitivity  of  the  release  altitude  and  other  model  pa¬ 
rameters  such  as  water  depth  and  water  temperature.  Altitude, 
when  varied  from  0  m  to  1000  m,  has  a  small  effect  on  burial 
depth  (relative  difference  of  18%).  When  a  more  realistic  upper 
limit  of  300  m  for  a  mine  laying  aircraft  is  applied,  the  relative 
difference  drops  to  9%.  Water  depth  affects  the  burial  depth  only 
if  a  mine  reaches  the  terminal  velocity  (in  this  case  about  20  m). 
Although  temperature  varies  the  water  density  up  to  3%  and  in 
turn  changes  the  viscosity  [29],  it  does  not  affect  the  burial  depth 
and  orientation. 

For  vertical  initial  falling  angle  (Fig.  6)  with  zero  rotation  rate, 
the  mine  is  heading  directly  downward,  resulting  in  the  maximum 
burial  depth.  When  the  release  height  is  150  m,  the  mine  burial 
depth  is  2.405  m  for  the  vertical  initial  falling  angle  and  0.359  m 
for  the  horizontal  initial  falling  angle.  Such  a  difference  in  burial 
depth  decreases  as  the  release  height  decreases.  When  the  release 
height  is  1.5  m,  the  mine  burial  depth  is  0.977  m  for  the  vertical 
initial  falling  angle  and  0.342  m  for  the  horizontal  initial  falling 
angle. 


Water  Depth  (m) 

0.5 1 - - - - - 

•§•0.45  ■ 

-C 

Cl 

&  0.4  ■ 

Is 

I  0.35  : _ : 

CD 

0.3’ - - - * - 

0  10  20  30 

(c)  Water  Temperature  (m) 

Fig.  5  Effect  of  (a)  release  attitude,  ( b )  water  depth,  and  (c) 
water  temperature  on  burial  depth.  Values  are  primarily  chosen 
to  represent  all  conditions  under  which  IMPACT25  and  IM- 
PACT28  may  be  used  (after  Chu  et  al.  [28]). 


5.3.2  Sensitivity  to  Sediment  Characteristics.  Figure  7  shows 
the  sensitivity  of  sediment  density  and  shear  strength  on  burial 
depth.  For  the  shear  strength  of  1  kPa  (extremely  soft  sediment), 
the  burial  depth  decreases  37%  for  sediment  density  varying  from 
1000  kg/m3  to  2000  kg/m3.  For  more  commonly  used  values  of 
shear  strength  (5-15  kPa),  the  sediment  density  has  very  little 
effect,  just  3.7%  (Fig.  7(a)).  As  shear  strength  increases,  the  in¬ 
fluence  of  sediment  density  reduces.  For  constant  sediment  den¬ 
sity  (1.5  kg  rrT3),  the  impact  burial  depth  drastically  reduces  from 
0.55  m  for  the  shear  strength  of  1  kPa  to  0.1  m  for  the  shear 
strength  of  10  kPa  (Fig.  7  (/?)). 

A  power  law  is  used 

S  =  apP  (42) 

to  represent  the  relationship  between  sediment  density  ( ps )  and 
shear  strength  ( S )  [30].  Here  the  two  coefficients  a  and  (3  are 
determined  experimentally.  Considering  a  homogeneous  sediment 
layer  (ps= const),  the  shear  strength  is  a  function  of  (a, (3).  For  a 
given  ps ,  increasing  a  or  /3  enhances  the  shear  strength  (S)  and  in 
turn  decreases  the  burial  depth.  These  features  are  well  simulated 
by  the  model  (Fig.  8). 

5.4  Strength  and  Weakness.  The  improvement  of  the  two- 
dimensional  model  versus  the  one-dimensional  model  is  its  capa¬ 
bility  to  predict  the  mine  rotation  in  the  (x,z)  plane,  i.e.,  ip2-  In  the 
two-dimensional  model,  the  momentum  equation  for  the  (x,z)  di¬ 
rections  and  the  moment-of-momentum  equation  for  the 
y-direction  are  used  to  predict  the  position  and  orientation  around 
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Fig.  6  Dependence  of  (a)  burial  depth  (m)  and  ( b )  height  protruding  (m)  on  release  altitude  (m)  and  water  depth  (m).  Height 
protruded  is  illustrated  here  to  clarify  the  levels  at  which  these  parameters  become  less  influential  in  the  two-dimensional 
models  (after  Chu  et  al.  [28]). 


the  y-axis  (i.e.,  in  the  (x,z)  plane).  The  basic  physics  for  the  spiral 
motion  described  in  Fig.  4  is  included  in  the  two-dimensional 
model  but  strictly  in  the  (x,z)  plane  [31,32]. 

Since  the  mine  movement  is  strictly  in  the  (jc,z)  plane,  it  is  very 
hard  to  include  the  motion  of  fluid  in  the  two-dimensional  model 
because  it  is  impossible  to  lay  a  mine  in  the  same  direction  of  the 
fluid  velocity.  In  the  littoral  zone,  the  water  velocity  is  not  negli¬ 


gible.  The  application  of  the  two-dimensional  model  for  the  op¬ 
erational  use  is  limited.  Besides,  the  drag  coefficients  for  the  axial 
and  cross  directions  have  similar  dependence  on  the  Reynolds 
number  Re  and  the  aspect  ratio  ( S )  (see  Eqs.  (17)  and  (18)).  How¬ 
ever,  the  drag  coefficients  in  the  axial  and  cross  directions  are 
independent  [33].  Besides,  a  mine  drop  experiment  shows  that 
IMPACT25/28  overpredicts  (five  to  ten  times  larger)  the  mine 


Fig.  7  Effect  of  sediment  (a)  density  and  ( b )  shear  strength  on  burial  depth.  Density  change  only 
impacts  the  predicted  burial  depth  in  very  soft  sediments.  As  expected,  shear  strength  has  a  dramatic 
impact  on  predicted  burial  depth  (after  Chu  et  al.  [31]) 
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alpha=5.4*l0 


-  rho=1000  kg/m3 

-  rho=1250  kg/m3 
rho=1 500  kg/m3 

-  rtio=1750  kg/m3 


Fig.  8  Effect  of  a  and  p  on  predicted  burial  depth  (m)  for  dif¬ 
ferent  values  of  sedimentary  density  (after  Chu  et  al.  [31]) 


burial  in  the  sandy  bottom  [9,31,32].  This  indicates  that  the  sedi¬ 
ment  dynamics  is  too  simple  in  the  two-dimensional  model. 


6  Three-Dimensional  Modeling  for  Cylindrical  Mines 

To  overcome  the  major  weaknesses  of  the  two-dimensional 
model,  i.e.,  (a)  environmental  fluid  assumed  motionless,  (b)  simi¬ 
lar  drag  coefficients  in  the  axial  and  cross  directions,  and  (c)  non- 
realistic  sediment  dynamics,  Chu  et  al.  [10-16,34-38]  modified 
IMPACT25/28  allowing  the  cylinder  to  move  in  three- 
dimensional  space.  These  changes  mandated  the  use  of  more  com¬ 
plicated  dynamical  system  than  IMPACT25/28.  The  three  coordi¬ 
nate  systems  (. F- ,  F-,  and  M-coordinates)  are  used  in  the  three- 
dimensional  modeling. 

6.1  Dynamical  System.  The  three  momentum  equations  (in 
the  ^-coordinate  system)  are  given  by 


u 

0 

V 

=  - 

0 

w 

_§  _ 

Ffr +  F/2 

Pn 


(43) 


and  the  three  moment-of-momentum  equations  (in  the 
M-coordinate  system)  are  written  in  vector  form 


Here,  (F^,M^)  are  the  buoyancy  force  and  torque,  and  (Fh,Mh) 
are  the  hydrodynamic  force  and  torque  including  the  drag,  lift,  and 
impact  (Fd,FhFi’,Md,MhMi).  The  vectors  (<o,Xl)  are  the  angu¬ 
lar  velocity  of  mine  and  M-coordinate  system, 

D  -  (45) 


The  first  term  in  the  right-hand  side  of  Eq.  (44)  is  an  apparent 
torque  (similar  to  the  Coriolis  term  in  earth  science)  due  to  the  use 
of  the  rotating  coordinate  system  (i.e.,  the  M-coordinate).  If  o)x 
=  0,  then  Q  =  o>.  The  apparent  torque  is  given  by 


-  2 j  •  (a  x  m) 

Jo  if  co1  -  0  (i.e., 11  =  of) 

^  —  2J 2  ^  i  W3J  m  +  2J 3  (o  j  if  oj  |  A  0 

(46) 


The  gravitational  force,  passing  the  COM,  does  not  induce  the 
moment.  In  the  M-coordinate  system,  the  moment  of  gyration 
tensor  for  the  axially  symmetric  cylinder  is  a  diagonal  matrix 


J  = 


Jx  0  0 

0  J2  0 

00/3 


(47) 


where  Jj,  J2 ,  and  J3  are  the  moments  of  inertia.  The  buoyancy 
force  induces  the  moment  in  the  j M  direction  if  the  COM  does  not 
coincide  with  the  COV  (i.e.,  ^A0), 

M/>=IFZ>I*C0S  i/yM  (48) 


6.2  IMPACT35-Cylindrical.  The  three-dimensional  model, 
usually  called  IMPACT35  (written  in  MATLAB),  is  the  third  gen¬ 
eration  of  the  Navy’s  mine  impact  burial  prediction  model.  For 
cylindrical  mines,  it  is  called  IMPACT3  5 -Cylindrical  or 
IMPACT35C.  It  was  developed  to  solve  the  six  scalar  equations 
(43)  and  (44)  to  obtain  (x,y,z, co1,co2, co3).  From  the  angular  ve¬ 
locity  ( (ol,(o2,(o3 ),  the  three  angles  determining  the  mine  orienta¬ 
tion  (i//j ,  i/f2,  ip3)  can  be  obtained, 


1  . 
dt 


dt 


■ 


d<h 

dt 


—  co3 


The  external  forcing  for  IMPACT35  consists  of  drag  force  Fd  and 
torque  Mj,  as  well  as  lift  force  F/  and  torque  M/.  All  these  exter¬ 
nal  forces  and  torques  are  calculated  in  the  F-coordinate  system. 
The  model  includes  cylindrical  and  noncylindrical  mines.  In  this 
subsection,  the  cylindrical  mine  is  taken  as  an  example  for  illus¬ 
tration. 

6.2.1  Drag  and  Lift  Forces.  Two  features  of  the  drag  and  lift 
coefficients  (Reynolds  number  dependent  and  asymmetry  in 
along-  and  across-mine  main  axis)  distinguish  IMPACT35  from 
IMPACT25/28.  The  drag  and  lift  forces  are  calculated  using  the 
F-coordinate  system.  Let  (Cdi,Cd2)  be  the  drag  coefficients 
along-  and  across-mine  main  axis  and  if\,f 2^/3)  be  the  added- 
mass  corrections  if 2^/3)  in  the  three  directions  of  the 
F-coordinate  system. 

The  total  drag  force  along  iF  (i.e.,  relative  flow  along  the  cyl¬ 
inder’s  main  axis)  is  calculated  by 

Fdl  =  Cldl(t)V  1,  |Vj(/)|  (49) 


J 


dco 

dt 


=  -  2J  •  (n  X  o>)  +  Mb  +  Mh 


(44)  Cdi  is  almost  independent  of  the  axial  Reynolds  number  (Re) 
when  Re>  104  but  dependent  on  the  cylinder’s  aspect  ratio  [33], 
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dx,  V2(x)  —  V2  —  (50) 


1.0 

0.75  +  #32.1934  +  0.09612/  S2 
1.15 


if  £>8 
if  8^  £>0.5 
if  0.5 


The  total  drag  force  along  jF  (i.e.,  relative  flow  across  the  cyl¬ 
inder)  is  calculated  by 


rL/2-x 

i  =  R  I  ( 

J  -L/2-x 


cjyj 


(i  +fi) 


where  V2  is  the  water-to-cylinder  velocity  at  the  surface  in  the  j F 
direction  and  an  empirical  formula  is  used  for  calculating  Cd2 
[39], 


J 


C/9  = 


1  .9276  +  8/Re 
1.261  +  16/Re 
0.855  +  89/Re 
0.84  +  0.00003  Re 
1.2-4  /S 
0.835 -0.35/£ 
0.7-0.08/£ 


if  Re  ^  12 
if  12  <  Re  ^  180 
if  180  <  Re  ^  2000 
if  2000  <  Re  «£  12,000 
if  12,000  <  Re  ^  150,000, 
if  12,000  <  Re  «£  150,000, 
if  12,000  <  Re  ^  150,000, 


£>  10 
2^  £<  10 
S<  2 


(51) 


1.875  -  0.0000045  Re  if  150,000  <  Re  «£  350,000 
^  1/(641550/Re  +  1.5)  if  Re  >  350,000 


The  drag  force  along  kF  is  calculated  by 

Frf3  = 


r»  Pw  F I  F I 


-L/2-X  r  0 

x2 dx  -  I  x2dx 
0  J  -L/2-X 


k  F 

(52) 


The  water-to-cylinder  velocity  determines  the  lift  force  [40], 


Ctl(t) 


rLU-x 

J  ^ 

J  -L/2-x 


(x)dx 


kF;  criW  -  \v2\ 


(53) 


c,= 


2(olR/V2  if  oi\RI\ 2  ’  4 

8  +  0.24(w1/?/y2  -  4)  if  <j>]R!V2  >  4 


(54) 


M„i  =  -  Cml  oj,iF,  Cml  =  TT/iLd1 


(55) 


Mj,= 


F I  F I 
—  0)2  1 0)2  | 


rL/2-x 

cd2R 

J  -L/2-x 


(1  +fr) 


(56a) 


where  fr  is  the  added-mass  factor  for  the  moment  of  drag  and  lift 
forces.  If  the  water-to-cylinder  velocity  or  the  cylinder  mass  dis¬ 
tribution  is  nonuniform  (^^0),  the  drag  force  causes  a  torque  on 
the  cylinder  in  the  kF  direction, 

rL/2-x 


l 


Cr]1R~ 


-L/2-x 


(1  +/r) 


(V2  -  oy2x)2xdx 


(56  b) 


The  lift  force  exerts  a  torque  on  the  cylinder  in  the  j F  direction, 

rL/2-x 


where  Q  is  the  lift  coefficient.  An  empirical  formula  is  used  for 
calculating  Q  [41], 


Ml2  = 


-l 


-L/2-x 


ctR—(y2  -  (4x)xdx 

fkr 


J  F 


(57) 


6.2.2  Drag  and  Lift  Torques.  For  an  axially  symmetric  cylin¬ 
der,  the  moment  of  the  hydrodynamic  force  in  \F  direction  is  not 
caused  by  the  drag  and  lift  forces,  but  by  the  viscous  fluid.  The 
moment  of  the  viscous  force  of  steady  flow  between  two  rotating 
cylinders  with  the  common  axis  is  calculated  by  [42] 

r  2  r  2 

M  =  4t7 — ^(ft>i  -  tuo) 


where  (rj  ,r0)  and  f oj  ,  ,&>0)  are  the  radii  and  the  angle  velocities  of 
the  inner  and  outer  cylinders;  yu  is  the  dynamic  viscosity.  The 
moment  of  the  viscous  force  on  one  rotating  cylinder  is  the  limit 
case  of  the  two  rotating  cylinders  as  r0^o°  and  co0=0.  The  mo¬ 
ment  of  the  viscous  force  around  iF  is  calculated  by 


Same  as  the  hydrodynamic  forces,  the  torques  along  the  j F  and 
kF  axes,  (Mdi,Md2,Mi),  are  calculated.  When  the  cylinder  rotates 
around  jF  with  the  angular  velocity  (o2,  the  drag  force  causes  a 
torque  on  the  cylinder  in  the  j F  direction, 


6.2.3  Interfacial  Treatment.  Computation  of  buoyancy  and  hy¬ 
drodynamic  forces  (F^Ffc)  and  torques  (M^,M/t)  is  more  compli¬ 
cated  for  a  cylinder  penetrating  through  air-water  and  water- 
sediment  interfaces  than  falling  through  a  single  medium  such  as 
water.  At  the  instance  when  the  cylinder  penetrates  into  an  inter¬ 
face,  three  situations  may  exist:  the  cross  section  is  a  complete 
ellipse  (Fig.  9(a)),  a  cutoff  ellipse  with  one  side  straight  line  (Fig. 
9(b)),  or  a  cutoff  ellipse  with  two  straight  lines  (Fig.  9(c)).  The 
interface  separates  the  cylinder  to  two  parts.  Each  part  contains  a 
noncylinder  D  and  a  subcylinder  C  (Fig.  10).  Let  (Lc,Ld), 
(HC,H^),  and  (Uc,Ud)  be  the  lengths,  surfaces,  and  volumes  of 
[ C,D ],  and  (hi,h2)  be  the  depths  of  the  two  sides  of  D  (Fig.  11). 
The  characteristics  of  the  geometric  parameters  (Lc,hi,h2)  are 
listed  in  Table  1.  The  COV  for  the  portion  [ C,D ]  is  called  the 
partial  COV  (PCOV).  With  this  treatment,  the  drag  and  lift  forces 
and  torques  at  the  two  parts  C  and  D  can  be  computed  separately 
[36]. 

6.2.4  Sediment  Resistance.  In  the  two-dimensional  models, 
IMPACT25/28,  the  pertinent  parameters  are  sediment  density  and 
bearing  strength.  The  bearing  strength  is  calculated  simply  by  ten 
times  the  shear  strength,  which  is  not  very  realistic.  In  fact,  when 
the  mine  impacts  and  penetrates  into  the  sediment,  it  will  create  a 
large  transient  pore  pressure  in  the  sediment  that  causes  ruptures 
in  the  sediment  and  influences  the  resistance  force  on  the  mine 
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Interface 


Interface 


(a)  (b)  (c) 

Fig.  9  Three  patterns  of  cylinder  penetration  with  the  cross  section  being  (a)  a  com¬ 
plete  ellipse  ( b )  cutoff  ellipse  with  one  side  straight  line,  and  (c)  a  cutoff  ellipse  with  two 
side  straight  lines  (after  Chu  and  Fan  [36]) 


surface  [11,36],  which  is  not  simply  ten  times  the  shear  strength. 
In  the  three-dimensional  model  (IMPACT35),  two  distinct  meth¬ 
ods  (delta  and  bearing  factor)  are  used  to  compute  total  sediment 
resistance  force  and  torque. 

(a)  Delta  method.  The  delta  method  is  developed  on  the  as¬ 
sumption  that  the  mine  pushes  the  sediment  and  leaves  space  in 
the  wake  as  it  impacts  and  penetrates  into  the  sediment.  This 
space  is  refilled  by  water  and  the  water  cavity  is  produced  (Fig. 


12).  At  the  instance  of  the  penetration,  the  total  sediment  shear 
resistant  force  on  the  mine  surface  is  calculated  by  [36] 


■=/  - 

a sed 


\i[lJ,G(V)S(z)]d(T 


-J  „[(aJ 


Ps(z')gdz' 


PwS^WS 


da  +  k^ps(z)\ 


gw  |  1  +eBdw\^ 
kp  ev  dt  / 


(58) 


where  in  the  right-hand  side  the  first  term  is  the  shear  resistance 
force,  the  second  term  is  the  buoyancy  force,  and  the  third  term  is 
the  pore  water  pressure.  S(z )  is  the  sediment  shear  strength,  G(V) 
is  the  impact  function,  V  is  the  mine  translation  speed,  ps(z)  is  the 
sediment  wet  density  (usually  obtained  from  the  sediment  data), 
(n,  r)  are  unit  vectors  normal  (outward  positive)  and  tangential  to 
the  mine  surface,  represents  the  vertical  coordinate  of  the 
water- sediment  interface  kp  is  the  permeability  coefficient 
(10-4  m  s"1,  [43]),  ev  (—0.50)  is  the  void  ratio,  and  B  is  the  length 
of  the  rupture  line.  The  step  function  A  is  defined  by 


[l,  v-n^0 

[o,  v-n^0 


(59) 


which  shows  that  the  sediment  buoyancy  and  shear  resistance 
forces  act  when  the  cylinder  moves  toward  it.  Let  \n  be  the  nor¬ 
mal  velocity.  The  tangential  velocity  is  represented  by 


Fig.  10  Illustration  of  PCOV  ( B~ ),  x1t  and  for  the  tail  part 
[C(1)5 0(1)]  for  the  case  in  Fig.  9(a)  (after  Chu  and  Fan  [36]) 


Fig.  11  Geometry  of  part  D(1)  (after  Chu  and  Fan  [36]) 


Table  1  Geometric  parameters  during  the  cylinder  penetration 
[36] 


4 

ft. 

h2 

Upper  and  lower  parts 

>0 

2  R 

0 

of  Fig.  9(a) 

Upper  part  of  Fig.  9(b) 

>0 

2  R 

0-2  R 

Lower  part  of  Fig.  9(b) 

0 

0-2  R 

0 

Upper  and  lower  parts 

0 

0-2  R 

0-2  R 

of  Fig.  9(c) 

\T=\-\n  (60) 

The  tangential  unit  vector  (r)  is  defined  by 

T=-pr  (61) 

H 

which  is  opposite  to  vT  (Fig.  13). 

The  sediment  resistance  torque  ( Ms )  is  calculated  by 


v 


Fig.  12  The  impact  (resistant)  force  exerted  on  the  part  of  the 
mine  surface  moving  toward  the  sediment  (after  Chu  and  Fan 
[36]) 
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Fig.  13  Momentum  and  angular  momentum  balance  for  mine 
penetration  through  the  water-sediment  interface 


Ms 


'^crsed 


(r  X  r)[A juG(V)S(z)]da 


I  (rXn)  (aJ 

•'"sed  L\ 


Ps(z')gdz’  +  pwg(zws-z ) 


77  ( gw  1  +  ev  dw 


en  dt 


+  (r  X  k)-ps(z)  —  + 


da 


(62) 


where  rpw  is  the  position  vector  (in  the  M-coordinate)  indicating 
the  location  of  the  cylinder’s  rupture  line. 

(b)  Bearing  factor  method.  The  bearing  factor  method  is  based 
on  the  fact  that  the  shear  resistance  force  (Fsr)  is  in  the  opposite 
direction  of  v  and  acts  on  the  mine.  Its  magnitude  is  proportional 
to  the  product  of  the  sediment  shear  strength  ( S )  and  the  rupture 
area  (A  is  projection  of  sediment-contacting  area  perpendicular  to 
the  velocity  V)  with  a  non-negative  bearing  factor  N  [44], 

Fsr  =  -  rN(p,v)SAda ,  N^O  (63) 

J  ased 


The  sediment  resistance  torque  includes  the  hydrodynamic  and 
shearing  resistance  torques, 


M  = 


{r  X  fh  -  [AT(/?,i;)SA]r  X  r}da 


(64) 


Here,  p  is  the  non-dimensional  penetration  depth  scaled  by  the 
diameter  (2 R).  The  sediment  density  and  shear  strength  S  in  Eqs. 
(63)  and  (64)  are  measured.  The  bearing  factor  increases  with  p 
and  decreases  with  the  decreasing  speed, 


N(p,v)  =  [jui1p^] 


1  +  X  logl 


(65) 


where  X  is  the  v -effect  parameter,  (piX  ,yu2)  are  the  p-effect  param¬ 
eters  [16],  and  ucri  is  the  critical  speed. 


6.2.5  Pseudocylinder  parametrization.  The  Navy  operational 
mines  are  usually  not  cylindrical.  It  is  important  to  develop  a 
model  with  more  general  shapes  such  as  with  nose  and  tail. 
Pseudocylinder  parametrization  was  proposed  for  noncylindrical 
mines  [35].  For  a  near-cylindrical  mine  with  nose  and  tail  falling 
through  a  single  medium  or  multiple  media,  the  buoyancy  force 
and  torque  are  relatively  easy  to  calculate.  But  the  hydrodynamic 


Fig.  14  Mine  with  nose,  tail,  and  cylindrical  body  (after  Chu 
and  Fan  [35]) 


forces  (lift,  drag)  and  torques  are  difficult  to  compute.  A  feasible 
way  is  to  transform  a  mine  with  nose  and  tail  to  a  cylindrical  mine 
(i.e.,  called  the  pseudocylinder  parameterization).  An  axially  sym¬ 
metric  mine  usually  consists  of  three  parts:  cylindrical  body  with 
radius  of  R,  nose,  and  tail  (Fig.  14).  The  lengths  of  the  mine,  nose, 
and  tail  are  L,  Ln,  and  Lt.  A  pseudocylinder  is  defined  with  the 
following  features:  same  radius  ( R )  of  the  mine’s  cylindrical  body 
and  the  same  volume  as  the  original  mine  (Fig.  15).  It  consists  of 
three  parts:  original  cylindrical  body,  and  equivalent  cylinders  for 
nose  and  tail.  Let  (n,nn,n?)  be  the  volumes  of  the  mine,  nose, 
and  tail.  The  equivalent  cylinder  has  length 


for  the  nose,  and 


JL 

7 tR2 


for  the  tail.  Let  (cc,cm)  be  the  mine’s  midpoint  on  the  main  axis 
and  the  COM  position,  and  let  cev  be  the  COV  of  the  pseudocy- 
lindrical  mine  (Fig.  15).  The  gravity  is  downward  and  passing 


Fig.  15  Location  of  cv,  cev,  and  cm.  Here,  e  is  the  distance 
between  cv  and  cm\  \  is  the  distance  between  cev  and  cm  (after 
Chu  and  Fan  [35]). 
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Table  2  Physical  parameters  of  the  model  mines  in  the  NSWC-Carderock  experiment  [47] 


Mine 

Mass  (kg) 

P 

(103  kg  m-3) 

L 

(m) 

7, 

(kg  m2) 

h  (j3) 

(kg  m2) 

X 

(m) 

1 

16.96 

1.60 

0.505 

0.0647 

0.356 

0 

2 

22.27 

2.10 

0.505 

0.0806 

0.477 

0 

3 

34.93 

1.60 

1.010 

0.1362 

2.900 

0 

4 

45.85 

2.10 

1.010 

0.1696 

3.820 

0 

5 

45.85 

2.10 

1.010 

0.1693 

3.940 

0.0045 

6 

45.85 

2.10 

1.010 

0.1692 

4.570 

-0.077 

through  cm.  The  buoyancy  force  is  upward  and  passing  through 
cev.  Let  Si  be  the  distances  between  cc  and  cm, 

Lfi  ~  Lne  Lt  —  Lte  ( 


Let  s2  be  the  displacement  from  cc  to  cm  that  is  easy  to  determine 
if  COM  is  given.  Let  x  be  the  displacement  from  cev  to  cm,  that  is, 
calculate 

X=si  +  s2  (67) 

Both  x  and  s2  can  be  positive  and  negative.  The  positive  values 
refer  to  nose-down  case,  i.e.,  the  point  cm  is  lower  than  the  point 
cev  for  positive  x  and  the  point  cc  is  lower  than  the  point  cev  for 
positive  s2. 

6.3  Mine  Impact  Burial  Experiment  for  Model 
Verification.  The  value  added  of  three-dimensional  model  (IM- 
PACT35)  versus  two-dimensional  model  (IMPACT25/28)  is  veri¬ 
fied  by  several  recent  mine  impact  burial  experiments  (Table  2): 
Mine  Impact  Burial  Experiment  (MIBEX)  at  Monterey  Bay  on 
May  22,  2000  [9,31,32],  Mine  Drop  Experiment  (MIDEX)  at  NPS 
swimming  pool  in  June  2001  [45,6,10,13,16,34],  MIDEX  at 
NSWC-Carderock  Explosion  Test  Pond  on  September  10-14, 
2001,  and  MIDEX  at  NSWC-Corpus  Christi  [46-48].  The 
MIDEX  at  the  Baltic  Sea  experiment  was  conducted  in  June  2003 
by  the  German  Federal  Armed  Forces  Underwater  Acoustic  and 
Marine  Geophysics  Research  Institute  [49,50]  with  the  full-size 
optical  mine,  which  is  allowed  to  free  fall  from  the  wench.  During 
these  experiments,  various  model  mines  (most  cylindrical)  were 
released  into  the  water.  The  mine  trajectories  were  recorded  by 
underwater  high-speed  video  cameras.  The  mine  burial  depths 
were  also  observed  by  the  diver  (in  MIB EX-Monterey  Bay)  and 
optical  instruments  (in  MIDEX-Baltic  Sea). 

The  two-dimensional  model  (IMPACT25/28)  and  three- 
dimensional  model  (IMPACT35)  are  integrated  using  the  same 
mine  parameters  (such  as  the  density  ratio,  length,  radius,  and 
distance  between  COV  and  COM)  and  mine  drop  initial  condi¬ 
tions  (speed  and  orientation)  as  in  the  mine  drop  experiments.  The 
value  added  of  IMPACT35  versus  IMPACT25/28  is  verified 
through  comparison  between  the  modeled  and  observed  mine  tra¬ 
jectories  and  burial  depths. 

6.4  Trajectory  in  Water  Column.  Detailed  verification  of 
trajectory  prediction  in  the  water  column  by  IMPACT35  and 
IMPACT25/28  has  been  reported  in  Ref.  [11]  using  data  from 
MIDEX  at  NPS  and  MIDEX  at  NSWC-Caderock.  For  near  hori¬ 
zontal  release,  the  3D  model  (IMPACT35)  simulated  trajectory 
agrees  well  with  the  observed  trajectory  with  the  same  travel  time 
(1.91  s)  of  mine  passing  through  the  water  column.  For  near 
45  deg  release,  the  3D  model  (IMPACT35)  simulated  trajectory 
and  travel  time  agree  well  with  the  observed  trajectory.  However, 
the  2D  model  (IMPACT28)  has  much  less  capability  to  predict  the 
cylinder  trajectory  in  the  water  column  with  near  horizontal  and 
45  deg  release.  For  near  vertical  release,  the  3D  model  (IM- 
PACT35)  simulated  trajectory  agrees  well  with  the  observed  tra¬ 
jectory  with  the  same  straight  pattern  and  the  same  travel  time 
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(1.83  s)  of  mine  passing  through  the  water  column.  However,  the 
existing  2D  model  (IMPACT28)  does  not  predict  the  travel  time 
well. 

6.5  Burial  Depth.  MIB  EX-Monterey  Bay  was  conducted  on 
the  R/V  John  Martin  on  May  23,  2000  [9,31].  The  barrel  with  a 
density  ratio  of  1.8  was  treated  as  model  mine  and  released  hori¬ 
zontally  while  touching  the  surface.  The  initial  conditions  are 

x(0)  =  y(0)  =  z(0)  =  0>  M(0)  =  v(0)  =  w(0)  =  0 

40)  =  90deg,  =  =  cof  =  <40)  =  <40)  =  0  (68) 

This  would  be  to  eliminate  any  chance  of  inertial  effects  caused 
by  uneven  introduction  into  the  air- sea  interface.  This  also  set  the 
initial  velocity  to  zero.  The  barrel  was  released  17  times.  The 
diver  would  snap  the  quick-release  shackle  on  the  barrel  and  then 
dive  down  to  conduct  measurements.  The  average  depth  of  the 
water  was  13  m.  Since  it  was  uncertain  the  path  the  barrel  would 
follow,  both  the  releasing  diver  and  a  second  safety  diver  would 
stay  on  the  surface  until  after  the  barrel  had  dropped.  Once  reach¬ 
ing  the  bottom,  one  diver  would  take  penetration  measurements 
using  a  meter  stick  marked  at  millimeter  increments  while  the 
other  would  take  a  gravity  core.  After  17  drops,  the  divers  began 
to  run  out  of  air  and  results  were  not  varying  greatly  so  the  deci¬ 
sion  was  made  to  end  the  experiment.  The  gravity  cores  were 
taken  immediately  to  the  USGS  Laboratories  in  Menlo  Park,  CA 
to  get  the  sediment  density  and  shear  strength  profiles  (Fig.  16). 

For  sediment  resistance  force,  the  two-dimensional  model 
(IMPACT25/28)  uses  ten  times  shear  strength  as  the  bearing 
strength.  The  three-dimensional  model  (IMPACT35)  uses  the 
delta  or  bearing  factor  method  (see  Sec.  6.2.4).  After  running  the 
two  models  (IMPACT35  and  IMPACT25/28)  for  each  gravity 
core  regime  [ps(z)  ,S(z)]  from  the  initial  conditions  (69),  the  burial 
depths  were  compared  with  measured  burial  depth  data  (Fig.  17). 


s, 

£ 
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Fig.  16  Mean  sediment  density  ps(z)  and  shear  strength  S(z) 
profiles  in  the  Monterey  Bay  collected  during  the  cylinder  drop 
experiment  on  May  31,  2000 
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Fig.  17  Observed  (MIBEX-NPS)  and  predicted  (IMPASCT25/28 
and  IMPASCT35  with  Delta  method)  burial  depths:  (a)  direct 
comparison  and  ( b )  scatter  diagram.  Note  that  the  two- 
dimensional  model  (IMPACT25/28)  predicts  the  burial  depth 
five  to  ten  times  larger  than  the  observed  depth,  and  IMPACT35 
with  Delta  method  performs  much  better  than  IMPACT/28  (after 
Chu  and  Fan  [36]). 


As  evident,  IMPACT35  improves  the  prediction  capability.  The 
2D  model  (IMPACT25/28)  overpredicts  the  actual  burial  depth  by 
an  order  of  magnitude  on  average.  However,  the  3D  model 
(IMPACT35)  predicts  the  burial  depth  reasonably  well  without 
evident  overprediction.  Since  the  gravity  cores  were  taken  for 
approximately  2-3  m  from  the  impact  location,  several  cores 
were  taken  for  each  drop.  This  allowed  an  average  to  be  calcu¬ 
lated  in  order  to  yield  more  accurate  data  for  each  drop.  Recently, 
Chu  and  Fan  [11]  compared  the  delta  and  bearing  factor  methods 
in  IMPACT35C  for  sediment  resistance  using  the  Baltic  Sea  ex¬ 
periment  data  and  found  that  the  bearing  factor  method  is  better 
than  the  delta  method. 

7  Modeling  of  Operational  Mines 

7.1  Shape  Effect.  Main  limitation  of  IMPACT35C  is  its  uti¬ 
lization  for  cylindrical  and  near-cylindrical  mines  only.  Shape  is  a 
significant  issue  if  the  model  is  used  operationally  because  the 
most  popular  mines  such  as  Rockan  and  Manta  are  far  from  cy¬ 
lindrical  (Fig.  18).  To  model  the  maneuvering  of  noncylindrical 
mines  in  water  column,  the  most  important  issue  is  to  determine 
the  hydrodynamic  (drag  and  lift)  force  and  torque  since  there  is  no 
existing  formula  for  calculating  the  drag  and  lift  forces  and 
torques  for  noncylindrical  objects.  Shape  effect  is  more  tenuous 


(a)  (b) 

Fig.  18  Operational  mines:  (a)  Manta,  and  ( b )  Rockan.  Here, 
the  Manta  is  an  anti-invasion  bottom  mine,  produced  primarily 
by  the  Italian  firm  Whitehead  Alenia.  It  is  shaped  as  a  frustum 
with  a  GRP  casing,  triggered  either  acoustically  or  magneti¬ 
cally.  The  Manta  has  a  shelf  life  of  30  years  and  will  operate  for 
17  months  after  activation.  The  Rockan  (made  in  Sweden)  has 
a  gliding  shape,  which  allows  mine  laying  over  a  wide  area 
while  covering  the  minimal  distance;  its  low-profile  stealth 
shape  makes  it  difficult  to  detect.  Its  casing  is  also  constructed 
of  GRP. 


Table  3  Physical  parameters  of  Rockan  mines 


Diameter 

Height 

Weight 

Charge 

Operating  depth 


0.980  m 
0.440  m 
220  kg 

130  kg  (HBX-3) 
3-100  m 


Table  4  Physical  parameters  of  Manta  mines 


Length 

Width 

Height 

Weight 

Charge 

Operating  depth 


1.015  m 
0.800  m 
0.385  m 
190  kg 

105  kg  (Cemtex) 
105  kg 


due  to  the  lack  of  test  data. 

MIDEX-II  at  NPS,  conducted  in  September  2005,  is  a  continu¬ 
ation  of  mine  drop  experiments  with  operational  mine  shapes 
[51,52].  In  that  experiment,  the  overall  shape  of  the  mine  was 
noncylindrical  mines.  In  addition  to  a  sphere  and  semihemispheri- 
cal  Gumdrop  shape,  two  shapes  were  specifically  chosen  to  rep¬ 
resent  real  world  bottom  mines:  the  Manta  and  Rockan. 

The  Manta  (Fig.  18(a))  is  an  anti-invasion  bottom  mine,  pro¬ 
duced  primarily  by  the  Italian  firm  Whitehead  Alenia.  It  is  shaped 
as  a  combined  frustum  and  flat  cylinder  with  a  glass  reinforced 
plastic  (GRP)  casing,  triggered  either  acoustically  or  magnetically. 
The  Manta  has  a  shelf  life  of  30  years  and  will  operate  for 
17  months  after  activation. 

The  BGM-100  Rockan  (Fig.  18(Z?))  is  an  acoustic  and  magnetic 
mine,  produced  by  Sweden.  It  has  a  gliding  shape,  which  allows 
mine  laying  over  a  wide  area  while  covering  the  minimal  distance. 
The  low-profile  stealth  shape  makes  it  difficult  to  detect.  Its  casing 
is  also  constructed  of  GRP.  Both  the  Manta  and  Rockan  are  made 
even  stealtheir  by  having  anechoic  coatings  and  being  made  of 


Table  5  Physical  and  geometrical  features  of  model  mines 
used  in  MIDEX-II 


Model 

Mass 

Density 

Dimensions 

Manta  bottom 

1145.0  g 

1.615  g/cm3 

D  (Bottom) 

15.0  cm 

influence  mine 

D  (Top) 

7.0  cm 

shape 

H 

6.2  cm 

Scale 

1/6 

Distance  from 

+0.373  cm 

COM  to  COY 

Z-axis 

Rockan  bottom 

813.0  g 

1.388  g/cm3 

L 

16.0  cm 

influence  mine 

W  (Back) 

7.8  cm 

shape 

W  (Front) 

13.3  cm 

H  (Back) 

6.3  cm 

H  (Front) 

3.8  cm 

Scale 

1/6 

Distance  from 

0  cm 

COM  to  COY 

Generic 

1692.0  g 

1.335  g/cm3 

D 

13.0  cm 

spherical 

Distance  from 

0  cm 

shape 

COM  to  COV 

Gumdrop 

2815.0  g 

1.722  g/cm3 

D 

14.9  cm 

hemispherical 

H 

13.3  cm 

shape 

Distance  from 

0  cm 

COM  to  COV 
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Mass  -  1 145.0  g 
Density  -  1.615  g/cm3 


Mass  -  813.0  g 
Density  -  1.388  g/cm3 


Fig.  19  Side  and  top  views  of  the  model  mines  with  different  shapes  for  MIDEX-II  experiment;  (a) 
Manta  and  ( b )  Rockan.  The  construction  of  these  model  mines  consisted  of  a  three-part  product 
process:  prototype  development,  mold  construction,  and  test  shape  casting,  and  finishing.  This  pro¬ 
cess  was  necessary  to  facilitate  more  efficient  experimentation  and  to  reduce  the  production  cost  of 
the  experimental  test  shapes. 


nonferrous  materials  to  reduce  the  mines’  acoustic  and  magnetic 
signatures.  Tables  3  and  4  show  the  geometric  and  physical  char¬ 
acteristics  of  Rockan  and  Manta  mines. 

7.2  MIDEX-II.  The  major  objective  of  the  MIDEX-II  was 
the  collection  of  trajectory  data  for  operational  mine  shapes.  Each 
shape  was  dropped  just  above  the  surface  of  the  water  and  filmed 
with  a  pair  of  high-speed  cameras  as  the  model  mines  fall  through 
the  water  column.  Each  trajectory  was  then  converted  to  an  array 
of  Cartesian  coordinates  (E-coordinate)  and  analyzed  with  soft¬ 
ware  specifically  designed  to  work  with  the  high-speed  cameras. 

Mine  shapes  are  selected  based  upon  current  and  future  opera¬ 
tional  relevance.  A  collection  of  four  minelike  polyester  resin  test 
shapes  were  used.  These  shapes  include  sphere,  semihemispheri- 
cal  “Gumdrop”  shape,  a  scale  model  of  the  Manta  bottom  mine, 
and  a  scale  model  of  the  Rockan  bottom  mine.  The  spherical  mine 
is  selected  to  serve  as  a  “calibration”  shape  because  its  symmetry 
and  equal  weight  distribution  are  about  its  three  axes.  The  Gum- 
drop  is  similar  in  shape  to  but  higher  in  density  than  the  Sphere. 
The  Gumdrop  is  selected  to  act  as  a  kind  of  “traditional”  shape  of 
bottom  sea  mines,  though  no  mine  was  specifically  represented. 
Table  5  shows  the  physical  and  geometrical  characteristics  of  the 
four  model  mines.  Figure  19  shows  the  geometric  features,  hori- 
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zontal,  and  lateral  views  of  Matha  and  Rockan  model  mines. 

MIDEX-II  was  conducted  at  the  Monterey  Bay  Aquarium  Re¬ 
search  Institute  (MBARI)  Unmanned  Underwater  Vehicle  Test 
Tank  (Fig.  20(a)).  Enclosed  inside  a  large  building,  this  10  X  15 
X  10  m3  tank  is  filled  with  “standard  sea  water.”  This  water  is 
maintained  by  an  ozone  filtration  system,  with  no  impurities  that 
save  the  remnants  of  dye  placed  into  the  tank  several  weeks  prior 
to  the  experiment.  The  faint  coloration  has  no  effect  on  the  shape 
trajectories,  but  it  did  add  some  difficulty  illuminating  the  tank. 
Hence  the  video  data  quality  is  somewhat  degraded.  A  sliding 
bridge,  on  which  the  slanted  board  is  mounted,  spanned  the  width 
of  the  tank.  Figure  20 (b)  describes  the  measurements  of  the  tank 
and  placement  of  the  drop  zone,  cameras,  and  lighting. 

The  tank  is  used  to  simulate  the  littoral  operating  environment 
with  the  scaled  depth  ranging  to  54.9  m  (180  ft);  however,  the 
nature  of  the  viewing  windows  only  allows  data  collection  to  a 
scaled  depth  of  roughly  18.29  m  (60  ft).  The  tank  has  no  current 
and  no  wind  blowing  over  its  surface.  To  aid  in  shape  recovery,  a 
9X  13  m2,  2  cm  netting  was  mounted  to  a  constructed  1.9  cm 
diameter  PVC  piping  horizontal  grid  and  placed  at  the  bottom  of 
the  tank  out  of  camera  view.  At  the  end  of  a  run,  the  net  that 
contained  the  shapes  was  retrieved  using  a  series  of  pulleys  placed 
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Fig.  20  MIDEX-II  setup:  (a)  MBARI  test  tank  facility  (structure  above  water  is  a  movable  bridge),  ( b )  top  view  of  the  two 
video  cameras,  (c)  view  from  underwater  viewing  window,  and  ( d)  calibration  test  cross.  Here,  the  MBARI  test  tank  (10 
X 15X10  m3)  was  filled  with  “standard  sea  water.”  This  water  was  maintained  by  an  ozone  filtration  system,  with  no 
impurities  that  save  the  remnants  of  dye  placed  into  the  tank  several  weeks  prior  to  the  experiment.  A  sliding  bridge,  on 
which  the  slanted  board  was  mounted,  spanned  the  width  of  the  tank  (see  (a)).  Eight  viewing  windows  (c)  were  1.83  m  (6  ft) 
below  the  surface  of  the  water.  The  two  viewing  windows  used  were  selected  because  of  the  unobstructed  and  near 
perpendicular  view  to  the  mine  drop  spot. 


at  the  four  corners  of  the  tank.  Two  large  dark  4.57X5.18  m2 
tarpaulins  were  placed  along  the  tank  walls  behind  the  camera 
views  to  assist  the  cameras  and  software  with  distinguishing  the 
falling  mine  shapes  from  the  tank  background.  Eight  viewing  win¬ 
dows  (Fig.  20(c))  are  1.83  m  (6  ft)  below  the  surface  of  the  water. 
The  two  viewing  windows  are  selected  because  of  the  unob¬ 
structed  and  near  perpendicular  view  to  the  drop  spot. 

All  the  data  were  collected  digitally  using  a  network  of  high¬ 
speed  and  standard  video  equipment  and  computers.  Data  above 
the  surface  are  collected  using  a  pair  of  standard  commercially 
available  digital  video  camera,  mounted  on  tripods,  and  located  at 
the  end  of  the  pool  directly  in  front  of  the  testing  zone.  Both  top 
cameras  operated  at  a  30  Hz  frame  rate.  The  data  camera  uses  a 
narrow  view  lens  zoomed  to  focus  on  the  area  directly  between 
the  slanted  board  and  the  water  surface,  and  is  toggled  on  and  off 
between  test  runs.  The  second  camera  used  a  wide-angle  lens  and 
was  employed  to  record  a  video  log  of  the  experiment.  This  de¬ 
vice  ran  continuously  throughout  the  experiment. 

Subsurface  video  data,  used  to  determine  the  trajectory  of  the 
falling  shapes,  are  collected  using  a  pair  of  high-speed  Photron 
FASTCAM  PCI  digital  cameras  (Fig.  20(d)).  These  cameras  are 
mounted  on  tripods  in  two  separate  windows,  at  an  angle  of 
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70  deg  in  relation  to  one  another  so  as  to  provide  two,  near  or¬ 
thogonal,  views  of  the  drop  zone.  After  mounting  and  calibration, 
each  camera  station  is  covered  with  black  plastic  to  block  out  any 
light  source  beside  that  which  comes  from  the  field  of  view.  The 
cameras  are  synchronized,  calibrated,  and  connected  by  a  cen¬ 
trally  located  laptop  computer  via  high-speed  data  cables.  During 
testing,  the  cameras  are  operated  using  the  Photron  FASTCAM 
Viewer  software  at  512X480  pixel  resolution  at  full  frame  and 
recording  rates  of  125  Hz.  To  facilitate  a  wider  field  of  view,  both 
cameras  are  fitted  with  wide-angle  lens.  All  data  are  recorded 
digitally  on  a  standalone  200  Gbytes  hard  drive  during  the  test 
phase.  Additionally,  to  enhance  the  quality  of  the  data  during  test¬ 
ing,  the  installed  tank  lighting  system  is  turned  to  its  maximum 
setting  and  a  pair  of  1000  W  high  intensity  photography  lights  are 
mounted  and  used  above  the  surface. 

The  mine  drop  experiment  consists  of  releasing  each  shape  ver¬ 
tically  from  about  0.3  m  (1  ft)  above  the  surface  of  the  water.  The 
entry  of  each  shape  into  the  water  is  recorded  by  the  two  above¬ 
surface  video  cameras.  All  subsurface  data  collection  is  facilitated 
by  the  two  FASTCAM  PCI  high-speed  cameras.  The  subsurface 
digital  data  are  analyzed  by  3D  motion  analysis  software  to  deter- 
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Fig.  21  Examples  of  high-speed  film  frames  for  model  mines. 
Here,  the  commercially  available  3D  motion  analysis  software, 
maxtraq,  was  the  primary  tool  utilized  to  perform  this  function. 


mine  the  trajectories  of  each  shape.  Overall,  we  have  15  drops  of 
the  Manta  mine,  14  drops  of  the  Rockan  mine,  9  drops  of  the 
Gumdrop,  and  13  drops  of  the  Sphere.  Initial  velocities  of  all 
shapes  as  they  entered  the  water  are  calculated  later  using  the 
maxtraq  motion  analysis  software. 

7.3  Data  Retrieval  and  Analysis.  Data  retrieval  has  been  ac¬ 
complished  following  all  the  experimental  test  runs  by  converting 
the  digital  video  imagery  (Fig.  21)  from  each  drop  into  an  array  of 
(x,y,z)  coordinate  data.  Commercially  available  3D  motion 
analysis  software,  maxtraq,  was  the  primary  tool  utilized  to  per¬ 
form  this  function.  Initially,  the  software  is  calibrated  into  the  3D 
coordinate  reference  system  utilizing  the  pairs  of  calibration  im¬ 
ages  obtained  in  the  initial  phase  of  the  experiment.  Following  the 
calibration,  both  camera  views  were  time  synced  and  analyzed  to 
determine  the  actual  position  of  the  shape  in  the  (x,y,z)  coordi¬ 
nate  field.  Frame  by  frame  analysis  was  performed  with  the  soft¬ 
ware  for  each  view  by  manually  identifying  and  inputting  one  or 
two  marker  points  associated  with  the  model  mine’s  position  and 
orientation.  For  the  Sphere  and  Gumdrop  shapes,  one  marker 
point  is  used  to  identify  the  lowest  position  of  the  shape.  For  the 
Manta  mine,  Point  1  (xi,yi,z\)  and  Point  2  (x2,y 2,z2)  are  chosen 
as  the  centers  of  the  bottom  diameter  and  top  diameter.  Examples 
of  temporally  varying  data  of  [xi(t)  ,yi(t)  ,z\(t)  ,x2(t)  ,y2(t)  ,Z2(t)] 
for  the  model  Manta  mine  are  listed  in  Table  6.  From  these  data, 
the  COM  location  (x,y,z)  and  the  orientation  (if/2,if/3)  are  calcu¬ 
lated  (Table  7). 

For  the  Rockan  mine,  Point  1  (xi,yi,zi)  is  selected  as  the  cen¬ 
ter  of  the  narrow  edge  at  the  “electronics  cylinder,”  and  Point  2 
(x2,y2,z2)  is  chosen  as  the  center  of  the  thicker  edge  of  the  shape 
opposite  to  Point  1 .  In  frames  where  a  marker  point  was  obscured, 
the  position  of  the  marker  (Point  1  or  Point  2)  is  estimated  visu¬ 
ally  based  on  the  previous  and  next  viewable  frame.  Following  the 
analysis  of  both  views,  the  automated  functions  of  the  software 
are  employed  to  compile  the  2D  images  into  a  calibrated  array  of 
3D  positional  data.  Examples  of  temporally  varying  data  of 
,yi(t)  ,Z\(t)  ,x2(t)  ,y2(t)  ,z2(t)]  for  the  model  Rockan  mine  are 
listed  in  Table  8.  From  these  data,  the  COM  location  (x,y,z)  and 
the  orientation  (i//2,if/3)  are  calculated  (Table  9). 

7.4  Trajectory  Patterns.  By  analyzing  the  2D  and  3D  plots 
of  each  shape  drop,  the  general  trajectory  patterns  of  the  four 
mines  were  obtained.  The  Sphere  and  Gumdrop  mines  are  hydro- 
dynamically  simple,  resulting  in  the  quickest  drop  times  in  the 
experiment.  The  Gumdrop  mine  has  the  fastest  mean  travel  time 
to  2.5  m  depth  (1.462  s).  The  Sphere  favors  a  straight-arc  trajec¬ 
tory  (0.62  probability)  over  a  curve-arc  trajectory  (0.38  probabil¬ 
ity).  The  Gumdrop  favors  a  curve  arc  (0.56  probability)  over  both 
the  straight  arc  (0.22  probability)  and  simple  slant  (0.22 
probability). 

The  Manta  and  Rockan  mines  have  more  complicated  shapes. 
Three  trajectory  patterns  (flat  spiral,  side  twist,  and  erratic)  are 


Table  6  Examples  of  data  for  model  Manta  mines  for  a  given  initial  velocity  (3.116  m/s) 


Time  (s) 

*1  (m) 

yi  (in) 

Zi  M 

*2  (m) 

yi  M 

Z2  M 

0 

0.2273 

0.6953 

-0.6371 

-0.2273 

-0.6953 

0.6371 

0.016 

-0.3519 

2.4886 

-2.821 

-1.1442 

-0.2139 

-3.5526 

0.032 

-0.3707 

0.7862 

-10.4696 

-0.7141 

-0.9996 

-9.4754 

0.048 

-0.8714 

0.531 

-16.3917 

-1.141 

-0.2135 

-13.5576 

0.064 

-0.5257 

0.1086 

-21.6666 

-0.8558 

-0.521 

-18.8631 

0.08 

-0.2535 

-0.3598 

-26.5455 

-0.9635 

0.2915 

-24.2111 

0.096 

0.0811 

-0.0842 

-30.5543 

-0.399 

-0.5669 

-28.8169 

0.112 

0.0192 

-0.8599 

-35.7727 

-0.8964 

-1.0305 

-33.9886 

0.128 

-0.3806 

-2.0507 

-39.9809 

-1.0781 

-1.6472 

-38.5195 

0.144 

-0.111 

-3.085 

-43.4884 

-0.236 

-2.4716 

-41.9996 
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Table  7  Center  of  mass  location  ( x,y,z )  and  orientation  (fc,  ^3)  for  model  Manta  mine  derived 
from  the  data  shown  in  Table  6 


Time  (s) 

x  (m) 

y  M 

z  (m) 

77/2-02 

fa 

0 

0 

0 

0 

-0.7165 

-1.8868 

0.016 

-0.748 

1.1373 

-3.1868 

0.2542 

-1.856 

0.032 

-0.5424 

-0.1067 

-9.9725 

-0.5003 

-1.7608 

0.048 

-1.0062 

0.1588 

-14.9747 

-1.2983 

-1.9182 

0.064 

-0.6907 

-0.2062 

-20.2648 

-1.3224 

-2.0538 

0.08 

-0.6085 

-0.0342 

-25.3783 

-1.1793 

2.3993 

0.096 

-0.1589 

-0.3256 

-29.6856 

-1.1973 

-2.3535 

0.112 

-0.4386 

-0.9452 

-34.8807 

-1.0897 

-2.9573 

0.128 

-0.7294 

-1.849 

-39.2502 

-1.0669 

2.6172 

0.144 

-0.1735 

-2.7783 

-42.744 

-1.1728 

1.7718 

identified  for  the  Manta  mine  (Fig.  22(a)).  Their  characteristics 
are  listed  in  Table  10.  For  example,  the  flat  spiral  pattern  shows 
the  Manta  mine  falling  with  its  bottom  side  basically  parallel  to 
the  X-Y  plane  and  following  a  spiraling  path.  Table  11  shows  the 
dependence  of  trajectory  pattern  and  travel  time  (to  2.5  m  depth) 
on  initial  speed  and  orientation  of  the  Manta  mine.  Three  trajec¬ 
tory  patterns  (flip-dive-flat,  flat- spin,  and  swoop-flat- spin)  are 


identified  for  the  Rockan  mine  (Fig.  22(b)).  Their  characteristics 
are  listed  in  Table  12.  For  example,  the  flip-dive-flat  pattern  shows 
the  Rockan  mine  flipping  once,  going  into  a  vertical  dive,  and 
settling  into  a  slowly  spinning  horizontal  orientation  for  the  re¬ 
mainder  of  the  drop.  Table  13  shows  the  dependence  of  trajectory 
pattern  and  travel  time  (to  2.5  m  depth)  on  initial  speed  and  ori¬ 
entation  of  the  Rockan  mine. 


Table  8  Examples  of  data  for  Rockan  mines  for  a  given  initial  velocity  (2.805  m/s) 


Time  (s) 

*1  (m) 

y\  M 

z  1  M 

*2  M 

yi  M 

Z2  M 

0 

-1.0736 

0.5296 

-5.166 

1.0736 

-0.5296 

5.166 

0.016 

-2.5339 

-0.2221 

-10.5596 

6.8934 

0.0791 

-0.5933 

0.032 

-4.0965 

-0.3357 

-12.7208 

6.3428 

-0.4706 

-5.2315 

0.048 

-6.2738 

-0.6902 

-14.0711 

6.3564 

-1.1554 

-10.9224 

0.064 

-7.6373 

-0.3879 

-13.5692 

6.2505 

-1.0931 

-15.9247 

0.08 

-8.1222 

-1.6365 

-15.0956 

2.8875 

-0.9055 

-19.5188 

0.096 

-7.7679 

-1.6508 

-15.1855 

2.2667 

-0.6023 

-21.3069 

0.112 

-6.2722 

-4.0571 

-22.3109 

0.1216 

3.0952 

-13.9761 

0.128 

-4.3282 

-5.2784 

-24.7814 

-1.2973 

2.809 

-17.453 

0.144 

-3.2107 

-5.5892 

-26.6419 

-2.2431 

3.5264 

-18.2012 

0.16 

-1.1347 

-5.8973 

-28.332 

-3.0412 

3.1642 

-19.8524 

0.176 

-2.7208 

-5.1987 

-30.0086 

-3.3042 

3.4589 

-23.1696 

0.192 

-2.4592 

-5.3611 

-32.0709 

-3.1931 

3.3691 

-24.8109 

0.208 

-2.4199 

-5.7808 

-34.2118 

-3.2777 

3.3184 

-26.3442 

0.224 

-1.7248 

-6.6254 

-36.7567 

-2.9388 

3.0883 

-27.4964 

0.24 

-1.8801 

-6.1255 

-36.9698 

-3.2004 

2.8093 

-30.25 

Table  9  Center  of  mass  location  ( x,y,z )  and  orientation  (faifa)  for  model  Rockan  mine  de¬ 
rived  from  the  data  shown  in  Table  8 


Time  (s) 

x  (m) 

y  M 

z  (m) 

77/2-02 

0 

0 

0 

0 

-1.3431 

-0.4583 

0.016 

2.1797 

-0.0715 

-5.5765 

-0.8129 

0.0319 

0.032 

1.1232 

-0.4032 

-8.9761 

-0.6223 

-0.0129 

0.048 

0.0413 

-0.9228 

-12.4968 

-0.2442 

-0.0368 

0.064 

-0.6934 

-0.7405 

-14.747 

0.1678 

-0.0507 

0.08 

-2.6173 

-1.271 

-17.3072 

0.3813 

0.0663 

0.096 

-2.7506 

-1.1265 

-18.2462 

0.5453 

0.1041 

0.112 

-3.0753 

-0.4809 

-18.1435 

-0.7153 

0.8413 

0.128 

-2.8128 

-1.2347 

-21.1172 

-0.7036 

1.2122 

0.144 

-2.7269 

-1.0314 

-22.4216 

-0.7442 

1.4651 

0.16 

-2.0879 

-1.3665 

-24.0922 

-0.7414 

1.7782 

0.176 

-3.0125 

-0.8699 

-26.5891 

-0.6675 

1.6381 

0.192 

-2.8262 

-0.996 

-28.4409 

-0.692 

1.6547 

0.208 

-2.8488 

-1.2312 

-30.278 

-0.7108 

1.6648 

0.224 

-2.3318 

-1.7686 

-32.1265 

-0.7576 

1.6951 

0.24 

-2.5402 

-1.6581 

-33.6099 

-0.6397 

1.7175 
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Fig.  22  Trajectory  patterns  of  model  mines  during  MIDEX-II:  (a)  Manta  and  ( b )  Rockan  mines 


The  Manta  mine  tends  to  either  fall  with  side-twist  (0.40  prob¬ 
ability)  or  with  flat-spiral  pattern  (0.40  probability),  and  occasion¬ 
ally  to  fall  with  erratic  pattern  (0.20  probability).  The  Rockan 
mine  tends  to  either  fall  with  swoop- flat- spin  (0.50  probability)  or 
flip-dive-flat  (0.36  probability)  pattern.  Both  trajectory  types  show 
the  potential  gliding.  For  two  drops  (0.14  probability),  the  Rockan 
mine  goes  directly  into  a  flat  spin  (Table  14). 

For  comparable  initial  speeds,  the  mean  travel  time  (to  2.5  m 
depth)  of  Manta  mine  (3.703  s)  is  over  twice  that  of  either  the 
Sphere  or  Gumdrop  mines.  The  mean  initial  speed  of  Manta  mine 
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Table  10  Trajectory  patterns  of  the  model  Manta  mine 


Flat  spiral  The  mine  falls  with  its  bottom  side  basically  parallel  to 
the  X-Y  plane  and  following  a  spiraling  path. 

Side  twist  The  mine  falls  with  its  bottom  side  perpendicular  to  the 
X-Y  plane  with  the  top  side  turning  about  the  bottom 
in  an  alternating  clockwise  and  anticlockwise  motion. 

Erratic  The  mine  falls  in  a  flipping  combination  of  side  and  flat 

spiraling  and  twisting. 
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Table  11  Distribution  of  trajectory  patterns  for  15  drops  of 
Manta  mine 


Table  14  Probability  of  trajectory-pattern  occurrence  for  all 
mine  shapes 


Drop 

Initial  speed 
(m/s) 

Initial  orientation 

Travel  time  (s)  to 
2.5  m  depth 

Trajectory 

pattern 

1 

1.819 

+X 

3.360 

Side  twist 

2 

1.573 

+X 

3.956 

Flat  spiral 

3 

2.160 

+x 

3.248 

Side  twist 

4 

1.531 

+X,  -Y 

4.800 

Flat  spiral 

5 

3.116 

1+ 

>< 

l 

4.040 

Flat  spiral 

6 

1.653 

— 

4.544 

Flat  spiral 

7 

1.881 

-Y 

5.040 

Flat  spiral 

8 

1.670 

+X 

4.656 

Flat  spiral 

9 

3.486 

-X,  +Y 

3.120 

Side  twist 

10 

2.265 

-X,  +Y 

2.940 

Side  twist 

11 

2.398 

+X,  +Y 

3.216 

Erratic 

12 

2.114 

-X,  +Y 

2.816 

Side  twist 

13 

3.143 

+Y 

2.928 

Side  twist 

14 

3.152 

+X,  +Y 

3.176 

Erratic 

15 

3.199 

-X,  +Y 

3.712 

Erratic 

Average 

2.344 

3.703 

Drop 

Initial  velocity 
(m/s) 

Initial 

orientation 

Travel  time  (s) 
to  2.5  m  depth 

Trajectory  pattern 

1 

2.341 

-X 

4.048 

Flip-dive-flat 

2 

2.805 

+X,  -Y 

3.984 

Flip-dive-flat 

3 

1.594 

-Y 

5.472 

Flat  spin 

4 

1.066 

+X,  -Y 

5.824 

Swoop-flat  spin 

5 

1.796 

+X,  -Y 

6.176 

Swoop-flat  spin 

6 

2.213 

+Y 

4.848 

Flat  spin 

7 

2.597 

-X,  -Y 

4.912 

Flip-dive-flat 

8 

2.652 

-Y 

4.880 

Flip-dive-flat 

9 

1.377 

+X,  +Y 

4.672 

Swoop-flat  spin 

10 

2.378 

+X,  +Y 

4.138 

Swoop-flat  spin 

11 

2.018 

+X,  +Y 

4.868 

Swoop-flat  spin 

12 

2.289 

+X,  +Y 

3.156 

Flip-dive-flat 

13 

1.872 

+X,  +Y 

4.740 

Swoop-flat  spin 

14 

2.362 

+Y 

3.912 

Swoop-flat  spin 

Average 

2.097 

4.688 

Shape  type 

Trajectory  type 

Probability 

Manta 

Flat  spiral 

0.40 

(15  total  drops) 

Side  twist 

0.40 

Erratic 

0.20 

Rockan 

Flip-dive-flat 

0.36 

(14  total  drops) 

Flat  spin 

0.14 

Swoop-flat  spin 

0.50 

(2.344  m/s)  is  around  11%  faster  than  the  Rockan  mine 
(2.097  m/s).  However,  the  mean  travel  time  (to  2.5  m  depth)  of 
Manta  mine  (3.703  s)  is  22%  shorter  than  the  Rockan  mine 
(4.688  s). 

7.5  Scaling  Effect.  Scaling  mines  down  from  the  actual 
mines  does  introduce  potential  for  altered  trajectory  based  on 
changing  the  drag  coefficient.  For  true  scaling,  the  Reynolds  num¬ 
ber  is  identical  for  the  full  scale  and  the  1/6  scale  model  (such  as 
in  MIDEX-II).  It  can  be  then  inferred  that  the  kinematic  viscosity 
is  reduced  by  the  scale  factor.  Since  the  mine  would  be  falling  in 
the  same  fluid  medium  (water),  kinematic  viscosity  remains  the 
same.  Thus  the  speed  of  the  falling  mine  determines  how  much 


Table  12  Trajectory  patterns  of  the  model  Rockan  mine 

Flip-dive-flat  The  mine  flips  once,  goes  into  a  vertical  dive,  and 
settles  into  a  slowly  spinning  horizontal  orientation 
for  the  remainder  of  the  drop. 

Flat  spin  The  mine  immediately  settles  into  a  slowly  spinning 

horizontal  orientation  and  remains  so  for  the  whole 
of  the  drop. 

Swoop-flat  spin  The  mine  makes  a  “U”  swooping  motion  upon  entering 
the  water,  after  which  it  settles  into  a  slowly  spinning 
horizontal  orientation  for  the  remainder  of  the  drop. 


Table  13  Trajectory  patterns  for  Rockan  mine 


drag  the  mine  experiences  to  affect  its  trajectory.  In  MIDEX-II, 
falling  speeds  of  noncylindrical  mines  were  primarily  in  the  lami¬ 
nar  flow  regime.  In  that  regime,  the  drag  coefficient  could  change 
with  different  Reynolds  numbers,  but  the  speeds  were  of  sufficient 
magnitude  that  the  drag  coefficient  remained  in  the  0.5-0. 6  range, 
reducing  the  drag  coefficient  error.  A  small  amount  of  the  speeds 
was  in  the  turbulent  flow  region,  where  the  drag  coefficient 
changes  very  little  for  different  Reynolds  numbers,  minimizing 
the  drag  coefficient  error  there  as  well.  So  the  scaled  models  pro¬ 
vide  a  good  representation  of  the  full  scale  mine.  The  Rockan 
scale  model  did  have  a  specific  trajectory  issue  independent  of  its 
Reynolds  number,  however.  Scaling  the  Rockan  shape  to  1/6  the 
actual  size  of  the  operational  mine  reduced  its  ability  to  “glide” 
through  the  water  as  easily  as  the  true  Rockan,  affecting  the  “re¬ 
alism”  of  the  Rockan  shape  trajectory. 

7.6  IMPACT35  for  Operational  Mine  Shapes.  Following 
the  procedure  in  model  development  of  IMPACT35C,  three  coor¬ 
dinate  systems  are  also  used  to  model  the  falling  operational 
mines  through  the  air,  water,  and  sediment  phases  (Fig.  23):  earth- 
fixed  coordinate  (E-coordinate),  main-axis  following  coordinate 
(M- coordinate),  and  force  following  coordinate  (E-coordinate) 
systems.  All  the  coordinate  systems  are  three  dimensional,  or¬ 
thogonal,  and  right-handed.  The  origin  of  M-  and  E-coordinates  is 
located  at  the  COM  of  each  operational  mine.  IMPACT35  for 
operational  mine  is  used  to  solve  the  six  scalar  equations  (43)  and 
(44)  to  obtain  (x,y,z,<u i,ft>2,co3).  From  the  angular  velocity 
(oq, (02,0)3),  the  three  angles  determining  the  mine  orientation 
(if/i,if/2,if/3)  can  be  obtained.  To  do  so,  we  need  first  to  derive 
formulas  for  calculating  the  hydrographic  (drag/lift)  forces  and 
torques  for  the  operational  mines.  Here,  we  take  IMPACT35- 
Manta  as  an  example  for  illustration. 

7.6.1  Geometric  Features  of  Manta  Mine.  Manta  mine  is  a 
combination  of  frustum  and  flat  cylinder  (Fig.  24).  The  frustum 
has  a  minimum  radius  (Ej)  and  a  maximum  radius  (E2)  with  a 
height  of  Lx.  The  flat  cylinder  has  a  radius  of  R2  with  a  height  of 
L2  (L2<L{).  Let 


R?  —  R^ 


V  — 


Ri 


(69) 


be  a  geometric  parameter  of  the  Manta  mine.  The  volume  of  the 
Manta  mine  is  calculated  by 


n  =  —{r\ + RiR2 + r22)  =  ttlr\\  1  +  rj+  — 


(70) 


The  distances  from  the  center  of  the  smaller  side  of  the  frustum 
(O)  to  the  COV  ( lv )  and  the  center  of  hydrodynamic  forces  (COF) 
(If)  are  calculated  by  (Fig.  24) 

rLl  '  L2 


j 

l  =  — 


ir[R(x)]2xdx  +  ttR2L2\  L 


.j]  + 

1  2 


LjQl7r[R(x)fdx  +  7 tR22L2 


L 

=  —  +  crn 
2 


(71) 
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Fig.  23  Three  coordinate  systems  for  Manta  and  Rockan  shapes 


f 

Lf~ 


2  R(x)xdx  +  2R2L2\L{  +  — 


$L^2R{x)dx  +  2R2L2 


L 

=  -  +  a, 


7 


(72) 


where 


Fig.  24  Location  of  cv,  cf,  and  cm.  Here,  \v  is  the  distance 
between  cv  and  cm\  \t  is  the  distance  between  cf  and  cm.  Here, 
R,  and  R2  are  the  small  and  large  radii  of  the  frustum. 


R(x)=R1  +  j-(R2-R1) 

is  the  varying  radius,  L=LX+L2  is  the  total  length  of  the  frustum, 
and  the  two  parameters  crv  and  ay  are  defined  by 


<rv  =  VL  L 

(6  +  6  77  +  2?;2)  +  6— (1  +  rj)1 
L 1 


1  +  - 


2+\v 


oy=  rjL 


6(2  +  rj)  +  12(1  +  rj) 


h 

U 


(73) 


The  COF  is  the  location  where  the  resultant  drag  force  exerts  with 
the  zero  moment.  Determination  of  the  COF  location  is  based  on 
the  assumption  that  the  drag  force  exerts  uniformly  on  the  Manta 
mine  surface.  The  distance  between  COM  and  COF  is  calculated 
by 


Xf=av-af 


(74) 


Here  Xf  is  used  to  replace  x  in  Eqs.  (56a),  (56 b),  and  (57)  for 
calculating  the  drag/lift  torques. 
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Experiment  MantalD:  3  time:  3.63  Impact35  time:  4.37  Experiment  MantalD:  4  time:  5.3  Impact35  time:  3.95 


Fig.  25  Model-data  comparison  of  Manta  mines  maneuvering  in  water  column 


7.6.2  Drag/Lift  Coefficients.  The  Manta  mine  has  large  aspect 
ratio  (S=2R2/L>  1).  Similar  to  the  cylindrical  mines,  the  drag 
force  is  decomposed  into  along  and  cross  axis  components.  The 
drag  force  along  the  main  axis  has  two  different  types  of  calcula¬ 
tion  depending  on  the  direction  of  the  flow.  When  the  water  flows 
from  R2  to  Ri,  the  Manta  mine  is  treated  as  a  disk  with  the  radius 
(. R )  of  R2.  When  the  water  flows  from  R{  to  R2 ,  the  Manta  mine  is 
treated  as  a  disk  with  the  radius  ( R )  of  ^Y1/(ttL),  where  II  is  the 
volume  of  the  Manta  mine.  Thus,  the  Reynolds  number  (Re)  is 
calculated  using  Re =2 RU/v. 

Dependence  of  drag  coefficient  along  the  main  axis  (Cd\)  on 
Reynolds  number  for  a  disk  has  been  published  as  a  curve  using 
the  experimental  data  by  Crowe  et  al.  [33].  From  this  curve,  a 
semiempirical  formula  is  derived  to  calculate  the  drag  coefficient 
along  the  main  axis, 

r  ( l05\ 

1.5  -  ,  Re  <105 

\  Re  / 


C,1  —  ^ 


1.5 +  0.3 


Re  -  105 
140-  105/’ 


105  <  Re  <  140 


,  Re  -  140 

1.8  -  (1.8  -  1.18)  — — —I,  140  <  Re  <  10 


104  -  140 


1.18, 


104  <  Re 


(75) 


On  the  other  hand  the  treatment  of  the  drag  force  across  the 
main  axis  and  the  lift  force  is  relatively  simpler  since  the  obstacle 
of  Manta  mine  to  the  flow  is  circular  (the  same  as  the  cylinder).  It 
is  reasonable  to  use  “equivalent  cylinder”  method,  i.e.,  to  compute 
the  drag  coefficient  across  the  main  axis  ( Cd2 )  and  lift  coefficient 
(Q)  using  Eqs.  (51)  and  (54)  with  the  same  volume  (II)  and 
length  (L)  of  the  Manta  mine  as  the  equivalent  cylinder. 


7.6.3  Moments  of  Inertia.  The  moment  of  gyration  tensor  in 
the  M- coordinate  system  for  the  Manta  mine  (i.e.,  a  combined 
flustum  and  flat  cylinder)  is  a  diagonal  matrix  with  the  diagonal 
components  (moments  of  inertia):  Jh  /2,  and  /3, 

rLi  ^  1 

/,=  -pR\{l+xlLx)Adx+-p'TrR\L2 

Jo  2  2 

TrpLlR\(5  +  l0r]+\0rjl  +  5rr’+  rft)  1  . 

= - +  -ottR?L2 

10  2r  2  2 


J2-  J3-  JA  +  JB 


(76) 


where  ( JA,JB )  are  the  inertias  of  the  frustum  and  flat  cylinder 
around  COM.  Let  the  distance  between  the  center  of  the  smaller 
side  of  the  frustum  (O)  to  the  COM  be  represented  by  lm  (Fig. 
24).  The  inertia  components  (JA,JB)  are  calculated  by 

rL\  ^4  rL  i 

JA=  I  - -(1  +x/Ll)4dx+  I  TTpR\{\  +x/Ll)2(x  -  lm)2dx 

J  0  4  J0 


=  P  1  1  (5  +  10 rj+  10 rf  +  5  rf*  +  rf)  +  ttpR\l\{-  -  —  +  ~rj 


20 

4  L 


3  U 


v~r  +  v\—  + 


3  'L 


L\2  .  rf-  rfl  lm\  rf(  lt 


2  \L 


—  —  +—  — 


3  \L 


(77) 


JB=^RtL2+^R22L32+  ttPR22L2 


(78) 


7.6.4  Model-Data  Comparison.  The  basic  equations  for 
IMPACT35-Manta  are  Eqs.  (43)  and  (44)  with  the  updated  com¬ 
ponents  (/]  ,J2,Jf)  (76)  in  the  gyration  tensor  J.  The  hydrody¬ 
namic  force  ( Fh )  and  torque  ( Mh )  are  calculated  using  Eqs.  (49), 
(50),  (52),  (53),  (55),  (56a),  and  (56 b)  using  the  updated  drag 
coefficients  ( Cdi,Cd2 )  and  lift  coefficient  C/  for  the  Manta  mine. 
In  the  model  development,  the  nonlinear  instability  and  model 
sensitivity  should  be  studied.  Within  the  correct  physics,  there  is  a 
possibility  of  chaotic  behavior  in  the  model.  The  chaotic  features 
will  be  handled  by  the  instability  and  predictability  analyses. 

IMPACT35-Manta  is  integrated  using  the  same  Manta  mine 
parameters  and  mine  drop  initial  conditions  (speed  and  orienta¬ 
tion)  as  in  the  mine  drop  experiments  (MIDEX-II,  see  Sec.  7.2). 
Performance  of  IMPACT35 -Manta  is  verified  through  comparison 
between  modeled  and  observed  mine  trajectories,  travel  times,  and 
positions.  The  3D  model  (IMPACT35-Manta)  simulated  trajectory 
agrees  well  with  the  observed  trajectory  (Fig.  25).  Both  show  the 
same  maneuvering  pattern  and  the  same  travel  time  for  the  Manta 
mine  passing  through  the  water  column.  The  modeled  and  ob¬ 
served  vertical  COM  positions  [zmod(0>£obsW]  are  well  compared 
for  all  the  15  drops  in  MIDEX-II  (Fig.  25).  Table  15  shows  the 
comparison  between  IMPACT35-Manta  and  MIDEX-II  (15  drops 
of  Manta  mine)  on  the  mine  maneuvering  depth  and  correspond¬ 
ing  travel  time.  The  mean  maneuvering  depth  predicted  by  the 
model  is  2.947  m,  which  is  very  close  to  the  observation 
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Table  15  Verification  of  IMPACT35-Manta  using  MIDEX-II  data 


Mean 

Mean 

(MIDEX-II) 

(Model) 

BIAS 

RMSE 

Travel  time  (s) 

4.425 

4.175 

-0.25 

0.9593 

Depth  (m) 

-2.942 

-2.947 

-0.005 

0.0057 

(2.942  m)  (Fig.  26).  The  travel  time  predicted  by  the  model  is 
4.175  s,  which  is  a  little  faster  than  the  observation  (4.425  s). 

8  Conclusions 

Advances  in  the  mine  impact  burial  prediction  in  the  past  three 
decades  are  reviewed  in  this  paper.  The  one-dimensional  model 
(IBPM)  was  developed  to  predict  the  vertical  position  of  the 
mine’s  COM.  The  model  provides  useful  information  such  as 
higher  falling  velocity  for  vertical  release  than  for  horizontal  re¬ 
lease,  strong  dependence  of  water  vertical  impact  velocity  on  the 


attitude  for  light  mines  (not  for  heavy  mines),  and  strong  depen¬ 
dence  of  the  mine’s  falling  velocity  in  water  column  on  attitude, 
wet  weight,  mine’s  length,  and  mine’s  radius.  For  the  same  mine, 
the  falling  velocity  in  the  water  has  a  minimum  in  horizontal 
orientation  and  a  maximum  in  vertical  orientation.  The  major 
weakness  is  the  constant  falling  angle  assumption  through  a  single 
fluid. 

The  two-dimensional  model  (IMPACT25/28)  was  developed  to 
overcome  the  major  weakness  of  the  one-dimensional  model  and 
to  predict  the  COM  position  in  the  (x,z)  plane  and  the  rotation 
around  the  y-axis.  The  major  weakness  is  the  difficulty  to  include 
the  fluid’s  motion  into  the  model  because  it  is  impossible  to  lay  a 
mine  in  the  same  direction  of  the  fluid  velocity.  In  littoral  zone, 
the  water  velocity  is  not  negligible.  The  application  of  the  two- 
dimensional  model  for  the  operational  use  is  limited.  Other  weak¬ 
nesses  include  (a)  a  similar  dependence  of  the  axial  and  cross  drag 
coefficients  on  the  Reynolds  number  Re  and  the  aspect  ratio,  and 
(b)  crude  parametrization  for  sediment  bearing  strength. 

The  three-dimensional  model  (IMPACT35)  has  been  developed 
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Fig.  26  Model-data  comparison  of  Manta  mines  maneuvering  in  water  column.  Here,  the  dashed 
curves  represent  model  results,  and  the  solid  curves  are  observations. 
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to  predict  the  COM  position  in  the  (x,y,z)  space  and  the  rotation 
around  the  x-,  y-,  and  z-axes.  IMPACT35  shows  great  improve¬ 
ments  versus  IMPACT25/28  using  the  recent  mine  drop  experi¬ 
mental  data.  Furthermore,  the  mine  shape  effect  is  also  included 
into  the  three-dimensional  model  (IMPACT35-Operational).  The 
most  important  features  of  the  Manta  mine  are  well  predicted 
when  the  modeled  data  are  compared  to  the  data  collected  during 
MIDEX-II  at  NPS  for  operational  mine  shapes  such  as  Manta  and 
Rockan.  The  drag  and  lift  laws  for  noncylindrical  mines  are  also 
derived. 
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Re 
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lift  coefficient 

translational  lift  coefficient  (kg  s-1) 
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added-mass  ratios  for  drag  and  lift 
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rotational  drag  force  (N) 
buoyancy  force  (N) 
drag  force  (N) 

drag  force  in  the  F-coordinate  (N) 
lift  force  (N) 

lift  force  in  the  F-coordinate  (N) 
unit  vectors  in  the  F-coordinate 
unit  vectors  in  the  F-coordinate 
unit  vectors  in  the  M-coordinate 
moments  of  gyration  (kg  m2) 
moments  of  gyration  for  cylindrical 
part-/  (kgm2) 
length  of  the  cylinder  (m) 
torque  due  to  the  buoyancy  force 
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torque  due  to  the  hydrodynamic  force 
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torques  due  to  the  drag  force  in  the 
M-coordinate  (kg  m  s-2) 
position  vector  (in  the  M-coordinate)  of 
point  on  the  cylinder’s  surface 
radius  of  the  cylinder 
Reynolds  number 
translation  velocity  (ms-1) 
water-to-cylinder  velocity  (m  s-1) 
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Xv 
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component  of  \r  along  the  cylinder 
(m  s-1) 

component  of  Vr  perpendicular  to  the 
cylinder  (ms-1) 
water  velocity  (ms-1) 
effective  cross  section  area  (m2) 
dynamic  viscosity  of  the  water  (m2  s-1) 
kinematic  viscosity  of  the  water  (m2  s-1) 
volume  of  the  cylinder  (m3) 
density  of  the  cylinder  (kg  m-3) 
density  of  the  water  (kg  m-3) 
distance  between  COM  and  COV  (m)  in 
cylindrical  mines 

distance  between  COM  and  COF  (m)  in 

noncylindrical  mines 

distance  between  COM  and  COV  (m)  in 

noncylindrical  mines 

angles  determining  the  cylinders’ 

orientation 

angular  velocity  (s-1) 

angular  velocity  components  in  the 

M-coordinate  (s-1) 

angular  velocity  components  in  the 

F-coordinate  (s-1) 
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